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1.  INTRODUCTION 


1.1  Solar  luminosity  variation 

The  appearance  and  disappearance  of  sunspots  during  the  11-year  solar  cycle 
has  long  been  suspected  of  causing  slight  fluctuations  of  the  Sun’s  luminous  output. 
Decades  of  pioneering  efforts  by  Abbott  and  his  coworkers  (1963)  from  several  sites, 
however,  failed  to  demonstrate  a  convincing  actual  detection  of  solar  variability. 
Groundbased  observations,  compromised  by  the  necessity  of  accounting  exactly  for 
the  variable  amounts  of  sunlight  absorbed  and  scattered  by  Earth’s  atmosphere, 
simply  were  not  capable  of  revealing  tiny  intrinsic  variations  in  the  ‘‘solar  constant.” 

In  1980,  extremely  precise  measurements  of  total  solar  irradiance  from  space 
finally  revealed  small  (<.0.2% )  dips  of  the  total  solar  output  accompanying  the 
disk  passage  of  large  sunspot  groups  (Willson  et  al.  1981).  Near  solar  maximum,  a 
strong  peak  in  the  power  spectrum  of  total  irradiance  at  the  solar  rotation  frequency 
is  caused  by  the  sunspot  flux  deficits  in  active  regions  (Willson  1982).  A  model 
of  solar  irradiance  variation  reconstructed  from  historic  sunspot  and  faculac  data 
predicted  that  owing  to  the  flux  blocking  effects  of  sunspots,  the  solar  output  should 
vary  inversely  with  solar  activity,  reaching  a  maximum  at  sunspot  minimum  and 
inversely  (Hoyt  arid  Eddy  1982). 

Contrary  to  this  simple  scenario  of  irradiance  variation,  continuing  measure¬ 
ments  since  1980  showed  that  the  solar  constant  decreased  by  almost  0.1%  from 
sunspot  maximum  to  sunspot  minimum  (Willson  et  al.  1986;  Willson  1987).  A  new 
model  of  the  radiation  balance  between  sunspot  flux  deficits  and  faculae  and  bright 
magnetic  network  flux  enhancements  was  proposed.  In  this  model,  excess  radiation 
from  faculae  and  network  components  compensates  the  flux  blockage  from  sunspots 
to  produce  the  observed  decline  of  total  irradiance  between  solar  maximum  and 
solar  minimum  (Foukal  and  Lean  1987,  1988;  Lean  and  Foukal  1988). 

With  the  existence  of  measurable  solar  variability  now  firmly  established,  the 
question  of  long-term  irradiance  changes  and  their  effects  on  the  solar- terrestrial  en¬ 
vironment  becomes  more  than  academic.  Even  though  the  sensitivity  of  terrestrial 
climate  to  solar  variability  is  currently  thought  to  be  much  smaller  than  anthro¬ 
pogenic  trace  gas  and  volcanic  effects,  the  baseline  of  measured  solar  variations— 
less  than  one  1 1-year  sunspot  cycle— is  still  far  too  short  to  permit  a  comprehensive 
assessment  of  the  causes  of  long-term  climatic  change  ( Dickinson  1987;  Ramanathan 
1988).  Thus  continuing  irradiance  monitoring  is  imperative. 

The  study  of  relationships  between  solar  irradiance,  solar  activity,  and  ter¬ 
restrial  weather  systems  is  also  becoming  important.  Although  past  efforts  were 
unconvincing  (e.g.,  Abbott  1963),  more  sophisticated  analysis  reveals  significant 
links  between  the  solar  cycle  and  some  large-scale  meteorological  phenomena  (Lab- 
itzke  1987;  Tinsley  198N).  Particularly  intriguing  is  the  situation  prevailing  during 
the  so-called  “Maunder  minimum”  from  about  1645  until  1715,  when  sunspots  were 
virtually  absent  (Eddy  1976).  An  accompanying  period  of  unusual  cold  in  northern 
Europe  was  so  extraordinary  that  it.  has  become  known  as  the  “Little  Ice  Age.”  Now 
that  an  apparent  connection  seems  to  exist  between  solar  activity  and  irradiance,  the 
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suggestion  of  a  causa]  relation  between  these  two  very  disparate  phenomena — one 
solar  and  one  climatological  -  seems  hard  to  avoid. 

It  is  in  this  context  that  the  study  of  solar  variability  has  undergone  a  re¬ 
naissance.  The  Sun  is  now  recognized  as  providing  a  varying  input  to  the  Earth 
radiation  budget,  and  some  component  of  this  variability  has  observable  terrestrial 
effects  on  solar-cycle  time  scales.  But  to  characterize  the  nature  and  significance 
of  solar  variations  by  observing  the  Sun  itself  will  require  decades  and  many  solar 
cycles.  For  a  more  timely  understanding  of  past  and  present  solar  variability,  we 
turn  to  the  stars,  especially  those  similar  to  the  Sun  in  age.  mass,  and  composition. 

1.2  Activity  cycles  in  Sun  like  stars 

Sun-like  chromospheric-activity  cycles  were  discovered  by  0.  C.  Wilson  (1978) 
in  the  course  of  a  decade-long  study  during  which  he  monitored  the  strengths  of 
the  emission  cores  of  the  II  and  K  lines  of  ionized  calcium  in  a  group  of  solar-type 
stars.  It  was  Wilson's  work  that  gave  rise  to  the  concept  of  “solar  analogs,”  stars 
that  truly  mimic  the  Sun  in  all  aspects  of  their  spectral  output  and  variability. 

A  basic  indicator  of  the  magnetic  activity  of  these  stars  is  provided  by  the 
variability  of  Wilson’s  index  of  composite  H-(-K  line  strength,  5,  defined  as  the  ratio 
of  the  emission  in  the  line  cores  to  the  nearby  continuum  level.  Stars  with  small, 
constant  S  indices  were  used  as  nightly  calibration  standards.  Old  main-sequence 
stars  with  the  weak  S  indices  signalling  low  chromospheric  activity  showed  activity 
cycles  remarkably  similar  to  the  Sun's,  some  comparable  in  length.  Continuing 
Wilson’s  program,  Vaughan  et  al.  (1981)  and  Baliunas  et  al.  (1983)  found  that 
young,  chromospherically  active  stars  displayed  rotational  modulation  of  the  HK 
flux,  attributable  to  the  stellar  analog  of  solar  plage,  thus  accounting  for  some  of 
the  random  fluctuations  observed  at  lower  time  resolution  by  Wilson. 

By  subtracting  the  photospheric  component  of  the  S  index,  Noyes  et  al.  (1984) 
derived  a  new  parameter,  the  “chromospheric  emission  ratio,”  that  quanti¬ 

tatively  estimates  the  actual  fraction  of  the  star’s  total  luminosity  originating  in 
the  chromosphere.  In  the  following  discussion,  we  shall  show  that  the  magnitude  of 
photometric  variability  is  closely  related  to  log  R'HK. 

1.3  Photometric  variability  of  Sun-like  stars 

Measurements  of  the  continuum  brightness  of  young,  chromospherically  active 
stars  in  the  Hyades  open  cluster  showed  that  nearly  all  of  these  stars  cooler  than 
spectral  type  F7  are  demonstrably  variable,  modulated  like  the  Sun — though  with 
much  larger  amplitudes — by  the  disk  transit  of  spotted  regions  (Lockwood  et  al. 
1984;  Radick  et  al.  1987).  Axial  rotation  periods  that  recur  from  year  to  year 
have  now  been  determined  for  18  stars.  Repeated  observations  in  successive  sea¬ 
sons  provided  evidence  of  long-term  changes  as  well:  during  seasons  when  the  rota¬ 
tional  modulation  of  brightness  was  most  apparent,  the  mean  brightness  was  reduced 
slightly  and  vice  versa.  A  simple  interpretation  is  that  the  presence  of  spots  reduces 
the  time-averaged  stellar  brightness — flux  deficits  associated  with  spots  are  not  com¬ 
pensated  by  excess  radiation  elsewhere  in  the  stellar  atmosphere,  as  is  evidently  the 
case  for  the  Sun. 


However,  in  contrast  to  the  clear  picture  that  has  emerged  regarding  the  re¬ 
lation  of  chromospheric  activity  and  photometric  variability  in  young  solar- type 
stars,  essentially  nothing  was  known  about  luminosity,  variations  of  common  main- 
sequence  stars.  In  fact,  some  of  the  Wilson  stars  exhibiting  Sun-like  chromospheric 
activity  cycles  are  among  those  chosen  decades  ago  and  long  accepted  as  photomet¬ 
ric  standard  stars  on  the  basis  of  their  presumed  stability.  Nearly  all  these  stars 
have  appeared  to  be  rock  steady  in  their  light  output.  A  photometric  program  car¬ 
ried  out  by  Jerzykiewicz  and  Serkowski  ( 19Gb)  from  1955  to  19G6  involving  repeated 
measurements  of  thirteen  mostly  main-sequence  P,  G,  and  K  stars  uncovered  no 
variability  at  the  level  of  about  1%. 

The  present  photometric  study  of  stars  similar  to  the  Sun  was  undertaken  in 
198  t  to  look  for  analog's  of  solar  variation  in  the  mostly  unexplored  and  techni¬ 
cally  difficult  regime  of  stellar  variability  below  1%.  Finding  evidence  of  luminosity 
variations  .n  a  group  of  similar  stars— especially  ones  with  known  Sun-like  activity 
cycles-  would  help  characterize  the  historic  variability  of  the  Sun  and  its  possible 
modes  of  variability  in  solar  cycles  to  come.  With  better  instrumentation  than 
Jerzykiewicz  and  Serkowski  had  available  25  years  ago,  we  hoped  to  measure  solar¬ 
like  stellar  variability  at  levels  well  below  0.5%  on  rotational  as  well  as  solar  cycle 
time  scales.  In  this  report,  vve  describe  the  results  of  the  first  four  observing  seasons 
of  this  program. 

2.  EXPERIMENTAL  DESIGN 

2.1  Selection  of  stars 

2.1.1  Program  stars 

Thirty  ix  main-sequence  and  subgiant  F,  G,  and  K  stars  were  selected  for 
photometric  monitoring,  all  but  two  from  Wilson’s  (1978)  survey  of  chromospheric 
activity  in  91  main-sequence  stars.  These  plus  others  suggested  by  D.  K.  Duncan 
(private  communication)  are  currently  also  monitored  spectroscopically  in  a  con¬ 
tinuation  of  the*  1IK  program  at  Mt.  Wilson.  Our  list  includes  HK  standards, 
young  active  stars,  and  Mjlar-age  stars  with  observed  cyclic  chromospheric  activity. 
The  program  stars  are  listed  in  Table  1,  organized  together  with  the  comparison 
stars  required  for  differential  photometry  into  30  trio  and  quartet  groups.  Program 
stars  are  denoted  by  tabulated  log  f?'HK  or  <  S  >  values.  Four  stars  (IID82885, 
HD1U710,  HD12013G.  1ID 13 1156)  were  previously  monitored  by  Jerzykiewicz  and 
Serkowski  (1966);  se\«  ral  are  urby  or  (.'BY'  photometric  standards,  11K  standards, 
or  MK  spectral  .-Tandards  and  are  so  designated  in  the  St ds  column  of  Table  1. 

In  Table  1,  the  "intra"  column  indicates  detected  variability  in  each  season, 
and  the  "inter”  column  indicates  variability  of  the  seasonal  mean  magnitudes,  as 
discussed  further  in  and  5. 
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TABLE  1.  Observing  List  and  Variability  Summary  (continued). 
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['ABLE  1.  Observing  List  and  Variability  Summary  (continued). 
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SUPPLEMENTARY  NOTES  -  TABLE  1 

The  information  in  these  notes  is  derived  in  large  pari  from  a  bibliographic 
search  kindly  conducted  for  us  by  Dr.  W.  II.  Warren,  Jr.,  National  Space  Science 
Data  Center,  using  the  SIMBAD  astronomical  hterature  database.  We  acknowledge 
with  gratitude  this  extremely  useful  service. 

In  general,  all  the  program  stars  and  many  of  the  comparison  stars  are  F2-K7 
dwarfs  of  apparent  magnitude  4.5  <  V  <  7  and  consequently  lie  close  to  Earth, 
many  within  10  pc.  Many  of  the  stars  are  thus  included  in  high  proper-motion 
and  parallax  catalogues  of  Gliese,  Luyten,  Woolley,  Lowell,  Yale,  etc.  Information 
about  the  dynamical  aspects  of  any  particular  star  should  be  sought  in  these  sources. 
Membership  of  several  stars  in  dynamical  groups  such  as  the  Hyades  are  noted  in 
the  table.  The  Mount  Wilson  UK  data  for  nearly  all  the  Wilson  stars  have  been 
analyzed  by  Vaughan  et  al.  (1981)  and  Baliunas  et  al.  (1983,  1985). 

IID1835  =  HR  88  =  9  Ceti  =  BE  Ceti.  This  Wilson  star,  1985  MK  G2V  standard, 
and  Hyades  group  member  has  been  frequently  studied  as  a  solar  analog.  However, 
because  it  is  younger  and  much  more  chromospherically  active  than  the  Sun,  it 
has  not  proven  to  be  a  perfect  match  (cf.  Hardorp  1978;  Cayrel  de  Strobel  et  al. 
1981;  Hardorp  1982).  Chugainov  (1980)  found  rotational  variability  (AV  =  0.032, 
AJ 3  =  0.037)  with  a  period  of  7.655  days,  using  HD2488  as  the  comparison  star 
(included  in  our  trio  with  HD1835).  Campbell  and  Cayrel  (1984)  found  molecular 
lines  characteristic  of  sunspot  umbrae  in  high-resolution  spectra  that  suggested  spot 
coverage  of  ~3%  of  the  surface. 

IID 10476  =  HR493  =  107  Piscium  =  NSV  600.  This  Wilson  star  is  a  photometric 
standard  on  the  UBV  (fundamental),  uu6j/,  and  DAO  systems.  The  suspected 
variability  is  due  to  Gutierrez-Moreno  et  al.  (1966),  who  claim  a  range  of  0.12  in  V, 
which  is  clearly  discredited  by  our  observations  (the  star  is  nevertheless  a  long-term 
microvariable).  The  star  is  a  suspected  spectroscopic  binary,  but  unresolved  via 
speckle  interferometry  (Hartkopfand  McAlister  1984). 

II I)  13421  =  I1R635  =  64  Ceti.  A  Wilson  star  that  has  received  relatively  little 
attention  apart  from  analysis  of  the  Mount  Wilson  HK  data.  It  is  included  as  a 
standard  in  Perry  et  al.  (1987).  MK  classification  suggests  this  is  a  subgiant. 

HD13611  =  HR649  =  Ceti  =  65  Ceti  =  NSV  749.  This  luminous  marginal 
barium  star  is  a  1983  MK  and  DDO  standard.  It  was  suspected  of  variability  by 
Argue  (1966),  who  suggested  it  was  an  eclipsing  binary  of  amplitude  V  =  0.03  from 
two  nights’  data  (!).  Our  data  also  show  variability,  but  of  a  different  kind;  short¬ 
term  variation  at  the  3%  level  is  ruled  out.  It  was  observed  by  Wilson  (1982)  during 
a  quick-look  survey  of  HK  emission  in  late-type  giants  and  found  to  be  inactive. 
The  star  has  a  diameter  determined  via  lunar  occultation  (cf.  White  and  Feierman 
1987).  It  is  a  spectroscopic  binary  with  several  published  orbits  (most  recently 
Griffin  1982),  which  is  resolved  via  speckle  interferometry  (McAlister  et  al.  1984); 
the  companion  was  detected  using  IUE  and  interpreted  as  a  white  dwarf  (common 
among  Ba  stars)  by  Bohm-Vitense  and  Johnson  (1985).  The  last  three  observations 
have  not  been  correlated  directly. 
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Supplementary  Notes  continued. 

II D 1 8256  =  HR869  =  />  Arietis  =  46  Arietis.  A  relatively  low  HK  flux  Wilson  star. 
The  reality  of  its  spectroscopic  binarity  has  been  disputed  (cf.  Morbey  and  Griffin 
1987). 

HD23140  =  SAO  39072.  Included  in  D.  K.  Duncan’s  1984  list  of  stars  then  actively 
observed  in  11 K  at  Mount  Wilson. 

HD25893  =  ADS  2995  =  V491  Persei.  A  Hyades  gro'  p  star  having  common  proper- 
motion  with  the  brighter  Wilson  star  IID25998  (see  t  jlow).  It  was  added  to  program 
because  of  probable  variability.  HK  emission  was  noticed  by  Bk.  'man  (1983)  on 
objective-prism  plates.  This  led  Boyd  et  al.  (1984)  to  check  it  for  variability  using 
an  automated  telescope;  they  found  a  V  amplitude  of  0.03  mag  and  a  period  of  7.37 
days.  Our  data,  though  not  as  extensive,  suggest  a  period  of  about  11  days  (more 
concordant  with  expectations  from  other  Hyades  stars  of  this  spectral  type),  and 
1986/7  season  amplitude  Ay  =  0.05. 

HD25998  =  1IR1278  =  50  Persei.  A  Wilson  star  and  Hyades  group  member  having 
common  proper-motion  with  the  binary  HD25893  (above).  Aside  from  analysis  of 
Mount  Wilson  HK  data,  no  other  relevant  study. 

I1D35296  =  HR1780  =  111  Tauri.  A  bright  Wilson  star  whose  chromospheric  ac¬ 
tivity  is  well-studied  in  many  lines  from  the  UV  to  IR.  Boyd  et  al.  (1984)  used  this 
as  a  comparison  star  for  the  long-period  variable  119=CE  Tauri,  but  did  not  notice 
the  variability  found  in  the  present  observations. 

HD39587  =  HR2047  =  y1  Orionis  =  54  Orionis.  Among  the  most  thoroughly  studied 
Wilson  stars,  and  among  the  brightest  on  our  observing  list.  Also  a  Stromgren 
photometric,  zero- polarization,  and  1985  MK  standard.  Its  angular  diameter  has 
been  measured  via  lunar  occultation  (cf.  White  and  Feierman  1987). 

HD75332  =  HR3499.  An  active  Wilson  star  without  much  astrophysical  study  apart 
from  Mount  Wilson  HK  analysis.  It  is  included  as  a  “highly  accurate”  photoelectric 
standard  in  Khaliullin  et  al.' s  (1985)  WBVR  system,  but  with  only  two  nights’ 
observation,  they  missed  the  2%  interannual  variation  we  see.  UBV  measures  have 
been  made  by  a  Canadian  amateur  (Kaitting  1985),  who  also  found  no  variations 
(his  rms  dispersion  was  0.02). 

HD76572  =  HR3563  =  61  Cancri  =  ADS  7107.  A  low-flux  Wilson  HK  standard. 
Classed  as  F6V  by  Cowley  and  Bidelinan  (1979),  but  moderate-dispersion  classifi¬ 
cation  by  Abt  (1986)  shows  it  to  have  slightly  weak  metallic  lines  (F4V/F2met/F4 
G-band).  Observed  by  Radick  et  al.  (1983)  in  1978/9,  who  found  no  variability  at 
the  4  mmag  level  in  y. 
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Supplementary  Notes  continued. 

111)81809  =  HR3750  =  Burnham  ‘2530.  'This  Wilson  star  is  often  considered  a  solar 
analog  (see,  however,  Neckel  1986).  Observed  by  Radick  (t  al.  (1983)  in  1978/9, 
who  found  no  variability  at  'lie  4  tntuag  level  in  y.  Duqueunov  and  Mayor  (1988) 
present  the  result  of  nine  years  of  CORAVEL  radial  velocity  measures,  which  are 
consistent  with  the  35-year  visual  binary  orbit  of  Baize  (1985);  the  implied  mass 
ratio  suggests  a  G0V+G9Y  pair. 

HD82635  =  HR3b00  =  10  Leonis  Minoris.  This  chromospherically  active  giant  was 
chosen  as  a  comparison  star  for  HD82885  =  11  LMi  (see.  next).  A  rotation  period  was 
found  by  Skiff  and  Lockwood  (1986).  The  star  is  a  standard  on  the  uvby  system. 
It  was  observed  by  Wilson  (1982)  in  a  quick-look  survey  of  11 K  emission  in  cool 
giants,  where  the  activity  was  about  double  that  of  the  Sun.  About  the  same  time, 
Middelkoop  (1982)  found  a  similar  S  value. 

H  1)82885  =  HI13815  =  11  Leonis  Minoris  =  ADS  7441.  An  active  and  well- 
studied  Wilson  star.  The  star  is  a  fundamental  UBV  standard,  uvby ,  DAO,  zero- 
polarization,  and  "Lowell  10- Year”  (Jerzykiewicz  and  Serkowski  1966)  standard. 
Radick  et  al.  (1983)  found  no  variations  in  1978/9  at  the  4  mmag  level.  Never¬ 
theless,  Skiff  and  Lockwood  (1986)  found  a  rotation  period  with  an  amplitude  over 
3%  in  1984/5.  The  spectral  type  has  been  disputed.  In  the  original  UBV  paper 
(Johnson  and  Morgan  1953),  Morgan  classed  it  as  G8IV-V;  later  the  YBS  listed 
it  as  G8III.  However,  the  large  parallax  and  proper-motion  indicate  that  the  star 
cannot  be  other  than  a  dwarf.  More  recently,  Abt  (1981)  has  classed  it  as  G8V.  The 
star  has  a  mag.  13  M-dwarf  companion  which  has  an  orbital  period  near  200  years 
(lleintz  1988). 

IID103095  =  HR4550  =  NSV  5374  =  CF  Ursae  Majoris  =  Groombridge  1830.  This 
thoroughly  studied  Wilson  star  is  a  famous  high-velocity  halo  subdwarf  8  parsecs 
from  the  Sun.  It  has  the  third-highest  known  proper-motion,  amounting  to  7.1 
arcsec-yr-1 .  The  star  is  used  as  a  UBV ,  Stromgren,  Washington,  Lowell  “10- Year,” 
and  radial-velocity  standard.  Flares  have  been  reported  twice  by  parallax  observers 
(see  lleintz  1984  for  summary),  but  these  almost  certainly  arose  from  instrumental 
effects;  the  variable  star  designations  were  thus  assigned  prematurely,  and  we  find 
only  long-term  changes  in  our  measures.  There  is  an  extensive  literature  of  detailed 
spectroscopy  and  spectrophotometry  examining  various  aspects  of  the  star’s  low 
metal  abundance  (roughly  ten  times  less  than  solar).  The  star  will  pass  near  another 
field  star  (V  —  13.90,  b  —  y  —  0.27,  based  on  two  nights  differential  photometry 
relative  to  HD  103095)  around  the  turn  of  the  century. 

111)111812  =  HR4883  =  31  Comae.  This  star  is  a  rapidly  rotating  G-giant  that 
appears  to  be  a  member  of  the  Coma  cluster,  but  lies  in  the  llertzsprung  gap.  It 
is  a  fundamental  UBV  standard,  Stromgren,  and  1985  MK  standard.  Well-studied 
spectroscopically  because  of  its  unusual  properties.  The  rapid  rotation  and  similarity 
to  FK  Comae-type  variables  led  Boyd  et  al.  (1984,  1985)  to  search  for  variability 
using  an  automatic  telescope;  none  was  found  at  the  1%  level  (Strassmeier  and  Hall 
1988);  our  measures  show  significant  seasonal  and  interannual  variations.  Simon 
(1986)  sought  rotational  modulation  in  IUE  spectra  without  success. 
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Supplementary  Notes  cont immd. 

HD  112989  =  HR-1924  -  .47  Comae  -  ADS  8731.  This  chemically  peculiar  K-giaut  is 
a  1985  MK  standard  ( G9li  1  ('ll  2Fe  1  (  a  1 )  and  has  been  much  studied  as  a  member 
of  the  class  of  metal-weak  stars.  !  he  visual  binary  companion  is  very  faint  (V  = 
13)  and  not  important  to  the  present  photometry.  Eggen  (1973)  included  the  star 
in  a  list  of  young,  disk  population  giants.  It  has  twice  been  classified  as  cooler  than 
the  standard:  KlIIIp  (weak  C  band)  by  Schild  (1973)  and  Kill-Ill  by  Abt  (1981). 
The  star  was  observed  photoelectricnlly  by  Boyd  ct  al.  (1981,  1985);  no  variability 
was  found  (Strassmeier  and  Hall  198s ). 

HD  1 11710  =  HR  1983  =  yf  Comae  =  13  Comae  Berenices.  This  thoroughly  studied 
Wilson  star  is  a  relatively  inactive  nearby  dwarf  and  1985  MK  standard. 

HD115383  =  HR5011  =  59  Virginis  =  e  Virginis.  A  fairly  active  Wilson  star  and 
Stromgren  standard.  No  directly  relevant  study  apart  from  HK  analysis.  Not  pre¬ 
viously  known  to  be  variable. 

HD115404  =  ADS  88-11  =  SAO  100191.  A  chromospherically  active  pair  involving 
K1  and  Ml  dwarfs.  Variability  was  first  reported  by  Radick  ct  al.  (1983),  who 
found  a  range  of  ~3 %  and  a  period  of  18  days,  consistent  with  Mount  Wilson  1IK 
measurements. 

HD  120 130  =  HR5185  =  4  Bootis  =  r  Bootis  =  NSV  644-4  =  ADS  9025.  This  Wilson 
star  is  a  U BV  and  Lowell  “10-year"  standard,  but  has  been  suspected  of  variability 
for  many  years.  The  chromospheric  activity  is  relatively  weak  as  measured  in  the 
HK  lines  and  in  He  I  lines  (see  WolfT  ct  al.  1986).  We  find  a  small,  gradual  change  in 
mean  brightness  after  four  seasons  unrelated  to  the  originally  suspected  6  Scuti-like 
variations,  which  do  not  exist  in  our  data.  The  faint  ( V  -  10.7)  companion  is  an 
M3  dwarf  (Joy  and  Abt  1974). 

HD124570  =  HR  5323  =  14  Bootis  =  NSV  6597.  This  Wilson  HK  standard  was 
observed  for  short-period  variability  by  Breger  (1969),  who  found  no  variations  at  the 
3  mmag  level  during  2.4  hours  of  monitoring.  The  spectroscopic  binary  companion 
is  undetectable  either  via  secondary  lines  (Gomez  and  Abt  1982)  or  via  speckle 
interferometry  (most  recently:  Blazit  el  al.  1987). 

HD129333  =  SAO  16453.  A  relatively  faint  star  on  the  UK  program  at  Mount  Wil¬ 
son  in  later  years.  The  most  active  star  in  our  photometric  program.  Well  observed 
in  the  Geneva  system,  but  not  suspected  of  variability  (see  Rufener  1981).  HK 
emission  analysis  only  by  Soderblom  (1985)  and  Soderblotn  and  Clements  (1987). 
Photometric  activity  quite  variable,  ranging  from  0.03  0.09  mag  in  y  over  one  sea¬ 
son,  and  changes  of  0.05  mag  or  so  over  several  years. 

HD131156  =  HR554  !  -  37  Bootis  =  £  Bootis  =  ADS  9413.  Another  active  double 
star  long  suspected  of  variability;  we  find  clear  evidence  of  rotation  despite  its  being 
a  Lowell  “10-year  standard"  (Jerzykiewicz  and  Serkowski  1966).  The  secondary 
star  is  2.23  magnitudes  fainter  in  V  on  average,  though  both  stars  are  likely  to  be 
variable. 
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Supplementary  Notes  continued. 

HD131511  =  HR5553  =  NSV  6847.  This  star  was  chosen  as  a  comparison  for 
HD131156,  but  as  it  is  a  member  of  the  Hyades  moving  group,  it  was  not  surpris¬ 
ing  to  find  it  variable;  our  data  provide  a  rotation  period  of  about  10  days,  near 
what  is  expected  for  Hyades  K2V  stars.  The  spectroscopic  binary  (P~125  days)  is 
unresolved  via  speckle  interferometry. 

HD143761  =  1IR5908  =  15  Coronae  Borealis  =  p  Coronae  Borealis.  This  low  HK- 
flux  high-velocity  halo  star  is  a  Stromgren,  radial- velocity,  and  1985  MK  standard. 

HD1496G1  =  HR6171  =  12  Ophiuchi  =  V2133  Ophiuchi.  An  active  Wilson  star  first 
discovered  to  be  variable  by  Dorren  and  Guinan  (1982)  after  null  results  published  by 
Blanco  ct  al.  (1979),  wdio  observed  the  star  for  four  seasons,  and  Rybka  (1979),  who 
analyzed  relatively  noisy  observations.  Radick  et  al.  (1983)  also  found  no  variability 
at  the  5  nunag  level. 

HD152391  =  SAOT21921.  An  active  Wilson  star,  but  photometric  variability  varies 
from  nil  to  obvious  between  seasons.  A  light  curve  was  first  published  by  Chugainov 
(1976).  Variability  was  also  observed  by  Dorren  and  Guinan  (1982),  who  derived  a 
rotation  period,  and  by  Radick  et  al.  (1983). 

HD157856  =  HR6489.  A  low  HK-flux  Wilson  star. 

HD158614  =  HR6516  =  ADS10598.  Another  inactive  Wilson  star.  Keenan  and 
Pitts  (1980)  included  it  their  1980  list  of  MK  standards  as  G9IV-V  1161,  however 
Abt  (1981)  classed  it  as  G5V. 

II D 1 60346.  Relatively  inactive  Wilson  star  and  1983  MK  standard.  Blanco  et  al. 
( 1979)  found  it  to  be  constant  at  the  1%  level. 

II D 1 6 1 239  =  I1R6608  =  84  Herculis.  The  sole  true  giant  on  the  Wilson  list  (relatively 
inactive),  and  a  1985  MK  standard  (G2IIIb). 

111)162211  =  11R6644  =  87  Herculis.  This  slightly  variable  comparison  star  has  tra¬ 
ditionally  been  used  as  a  comparison  star  for  the  extensively  observed  UU-Herculis 
variable  89  Her;  the  variability  we  find  has  not  been  previously  suspected  (e.g., 
Fernie,  1986,  where  Stromgren  photometry  is  given). 

HD176095  =  HR7163.  A  relatively  hot  Wilson  star. 

HD182572  =  HR7373  =  b  Aquilae  =  31  Aquilae  =  NSV  11994.  This  Wilson  star  has 
attracted  attention  largely  because  of  its  high  velocity  and  its  strong-lined  spectrum. 
It  is  a  radial-velocity  (extensive  high-precision  observations  indicate  it  is  single)  and 
1985  MK  standard. 

HD185018  =  HR745G  =  ADS12670.  This  supergiant  was  selected  as  a  comparison 
star  before  we  were  aware  of  its  luminous  nature.  It  lies  just  outside  the  Cepheid 
instability  strip,  and  variability  has  not  been  demonstrated  before  the  present  ob¬ 
servations  (Eggen,  1985,  found  it  nonvariable  at  the  3  mmag  level  on  four  nights) 
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Supplementary  Notes  continued. 

HD185144  =  IIR7462  -  01  Draconis  =  a  Draconis  ~  NSV  12170.  This  bright, 
nearby  K  dwarf  is  a  UBV.  Stromgren,  and  1985  MK  standard.  The  suspected 
variability  (AT  =  0.07)  was  found  visually  and  can  be  discounted  from  our  data. 

1ID1900U7  =  SAO  125379.  A  cool,  active  Wilson  star.  A  marginal  detection  of 
variability  was  found  by  Dorren  and  Guinan  (1982),  but  there  were  insufficient  data 
to  make  a  firm  statement.  Radick  ct  al.  (1983)  also  came  to  an  uncertain  conclusion 
based  on  Stromgren  data  whose  y-filter  rms  dispersion  was  0.009  mag. 

111)201091  =  HR8085  =  61  Gygni  A  --  VT803  Cygni  =  ADS14636A.  This  famous, 
exhaustively  studied  star  (and  UBV,  Stromgren,  and  1985  MK  standard)  was  first 
clearly  observed  to  be  variable  by  Dorren  and  Guinan  (1982).  They  found  an 
amplitude  of  0.03  mag  in  the  blue  that  correlated  well  with  the  HK-derived  rotation 
period.  They  observed  near  the  cyclic  HK-finx  maximum  about  1981.  In  contrast, 
Perry  et  al.  ( 1987)  find  no  variability  on  18  nights;  likewise  our  observations  do  not 
show  this  short-period  variability. 

IID201092  =  HR80S6  =  61  Cygni  13  =  NSV  13546  =  ADS14636B.  A  standard  on 
many  systems:  UBV.  Stromgren,  DAO,  1985  MK,  and  others.  This  star  has  also 
been  suspected  of  variability  (AT  =  0.07),  which  is  ruled  out  by  the  present  obser¬ 
vations.  Neff  (1968)  claimed  "rapid  irregular  light  variations  or.  several  occasions.” 
Blanco  f  t  al.  ( 1979)  saw  no  variations  over  five  seasons;  Dorren  and  Guinan  (1982) 
had  insufficient  data  to  come  to  a  conclusion  about  variability;  finally,  Perry  ct  al. 

(  1987)  found  no  variability  over  18  nights. 

HD215427  =  SAO  52320.  A  new  semiregular  variable  among  our  comparison  stars 
for  HD2157Q4.  No  citations  in  the  SIMBAD  database. 

HD215704  =  SAO  52353.  A  relatively  faint,  low  HK-fiux  star  from  the  1984  Duncan 
observing  list. 

II D2 16385  =  HR8697  --  49  Pcgasi  =  a  Pegasi.  UBV  and  Wilson  IiK  standard. 
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Figure  1  shows  the  relationship  between  the  chromospheric  emission  ratio  (log 
R'hk)  and  B  —  V  color  for  the  program  stars.  The  designation  “Y”  or  “0”  refers 
to  the  classification  of  the  stars  as  “young”  or  “old,”  according  to  their  position  in 
an  [<  S  >,  B  -  lr]  diagram  (Vaughan  1980).  On  this  diagram  the  Sun’s  position 
is  bracketed  in  color  and  activity  level  by  program  stars  that  cover  the  full  range  of 
both  quantities.  In  these  parameters,  the  closest  solar  analog  is  HD81809,  a  G2V 
spectroscopic  binary.  The  location  of  the  Sun  on  Figure  1  suggests  that  it  may  be 
among  the  hottest  of  the  old  stars  having  a  measurable  activity  cycle.  At  the  lower 
left,  the  six  stars  denoted  by  filled  symbols  are  HK  standards  with  low  and  constant 
<  S  >  values;  two  of  them  (HD13421,  GOV;  HD143761,  GO+Va)  happen  to  be 
uvby  photometric  standards  as  well,  and  all  have  subsolar  activity  levels  according 
to  their  low  values  of  log 


Figure  1.  Chromospheric  emission  ratio  log  R'HK  as  a  function  of  B  —  V  color 
for  the  program  stars.  (Y)  -  young,  active  stars,  (O)  -  old  stars,  ( filled  circles) 
inactive  HK  standards.  The  designations  “old”  and  “young”  are  from  Vaughan 
(1980). 


2.1.2  Comparison  stars 

Optimal  differential  photometry  requires  the  measurement  of  two  or  more  field 
comparison  stars  of  comparable  apparent  magnitude  and  color.  A  minimal  configu¬ 
ration  is  thus  a  trio  containing  one  program  star  and  two  comparison  stars  located 
within  a  few  degrees  of  one  another  on  the  sky.  In  six  cases  initially,  quartet  groups 
were  formed  containing  two  program  stars  and  two  comparison  stars.  In  Table  1,  the 
entries  are  organized  into  the  resulting  30  such  groups.  When  more  than  one  star 


-  14  - 


Declination  (degrees) 


in  a  trio  exhibited  variability,  an  additional  comparison  star  was  added,  eventually 
making  several  new  quartet  groups. 

Figure  2  shows  the  location  of  the  groups  on  the  sky.  The  star  symbols  are 
program  stars.  In  most  cases,  the  group  members  are  so  close  together  that  their 
airmasses  are  very  nearly  the  same,  thus  minimizing  differential  extinction  effects. 
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Figure  2.  Location  of  the  trio  and  quartet  groups  on  the  sky.  (stars)  -  program 
stars.  Differential  airmasses  may  be  estimated  by  the  separate  scale  showing  the 
airmass  at  meridian  transit. 


2.1.3  Distribution  of  color  and  apparent  magnitude 

Figure  3  shows  that  the  distributions  of  apparent  V  magnitudes  and  B  —  V 
colors  for  the  program  stars  and  comparison  stars  are  similar,  though  the  comparison 
stars  are  on  average  about  0.5  mag  fainter  and  0.0 1  mag  redder  in  B  —  V . 

2.2  The  observational  facility 

All  observations  were  made  using  the  0.5-m  reflecting  telescope  of  Lowell  Ob¬ 
servatory,  located  on  a  wooded  inesa  at  an  elevation  of  2200  m  about  1  km  west 
and  100  m  above  downtown  Flagstaff.  The  telescope  is  a  manually  slewed  f/ 16 
Cassegrainian  housed  in  a  building  with  a  roll-off  roof;  it  was  originally  used  in  the 
1950s  by  H.  L.  Johnson  in  the  early  days  of  UBV  photometry.  Recent  photographs 
of  the  telescope  may  be  found  in  Lockwood  (1983)  and  Lockwood  and  Skiff  (1987). 


V  magnitude 


Figure  3.  Apparent  visual  magnitudes  and  measured  or  estimated  B  -  V  colors 
for  the  program  stars  {stars)  and  comparison  stars  (open  circles). 


Despite  the  proximity  of  city  lights,  the  sky  brightness  is  only  a  couple  of 
stellar  magnitudes  brighter  at  the  zenith  than  the  natural  celestial  background  and  is 
therefore  essentially  negligible  for  bright-star  photometry.  Local  aerosols  are  usually 
trapped  below  the  level  of  the  Observatory  by  a  nighttime  temperature  inversion,  a 
significant  factor  in  winter  when  woodstoves  are  in  use. 

Since  1974.  the  telescope  has  been  equipped  with  a  photoelectric  photometer 
of  conventional  design,  pulse-counting  electronics,  and  a  data-acquisition  microcom¬ 
puter,  currently  a  Digital  Equipment  Corporation  LSI-11/03  and  associated  stan¬ 
dard  interface  and  control  circuit  cards.  Data  are  displayed  on  a  video  terminal  and 
recorded  on  DEC  RX0‘2  flexible  disks  for  subsequent  analysis  on  the  Observatory’s 
VAX- 750  computer.  Recently,  funds  provided  by  this  contract  were  used  to  equip 
the  telescope  with  a  computerized  coordinate  display  system. 

The  photometer  utilizes  a  thermoelectrically  cooled  EMI  6256S  photomulti¬ 
plier  and  standard  intermediate-band  Strbmgren  uvby  interference  filters.  An  Sr90 
Cerenkov  source  provides  a  stable  light  output  used  to  monitor  the  nominal  perfor¬ 
mance  of  the  detector  and  electronics.  The  same  photomultiplier  and  filters  have 
been  in  constant  use  since  1971,  a  factor  contributing  significantly  to  the  long-term 
stability  and  precision  of  results  obtained  with  this  facility. 
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2.3  T ho  basic  observation  sequence 

Following  the  methodology  reconunetuled  to  us  by  Jerzykiewicz  15  years  ago, 
the  basic  data  unit  in  precision  differential  photometry  is  a  “cycle”  comprising 
sequential  observations  of  the  three  stars  of  a  trio  or  the  four  stars  of  a  quartet 
through  a  single  filter.  In  this  program  we  utilize  only  the  b  (172  nm)  and  y  (551 
nm)  filters  of  the  Stromgren  four-color  system.  Typically,  a  standard  nightly  obser¬ 
vation  consists  of  four  cycles  in  the  filter  order  y,b,b.y.  Within  eacli  cycle,  each  star 
is  measured  for  a  total  of  six  10-second  intervals  as  follows: 
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**  -  quartet  only. 

Subsequent  cycles  in  the  same  filter  begin  at  f;  subsequent  cycles  in  a  different 
filter  begin  at  * 

A  sample  complete  observation  is  given  in  the  Appendix.  Each  set  of  n  cycles 
is  preceded  and  followed  bv  a  dark  count  and  standard  light  source  sequence  that 
servos  as  an  archived  indicator  of  system  performance. 

An  important  element  of  the  experimental  design  is  the  equal  statistical  weight 
(  in  terms  of  total  integration)  of  each  stellar  measurement,  regardless  of  whether  it 
is  a  program  or  a  comparison  star.  The  various  pairwise  differential  combinations 
(i.e.,  star  1-star  2,  star  1-star  3,  ....  star  3-star  4)  are  thus  statistically  equiva¬ 
lent  except  for  time  order  and  the  splitting  of  the  observations  of  star  2  (normally 
a  comparison  star)  into  two  half- weight  sequences  beginning  and  ending  each  cy¬ 
cle.  Every  nightly  differential  magnitude  set  for  a  group  contains  an  internal  noise 
measure  provided  by  two  presumably  quiescent  comparison  stars.  The  most  stable 
pairwise  differential  magnitude  defines  the  night-to-night  noise  baseline,  which  in 
t  he  absence  of  variability  i-,  typically  in  the  0.002  0.003  mag  range. 

2.4  Expected  precision  levels 

Attaining  the  higher  possible  long-term  precision  is  a  major  goal  of  this  pro¬ 
gram.  If  they  mimic  the  Sun  exactly,  our  solar  analog  program  stars  will  have  rota¬ 
tional  amplitudes  under  0.5%  and  cycle  amplitudes  under  0.1%.  This  is  a  formidable 
challenge  to  conventional  stellar  photometry,  one  which  has  few  precedents  in  the 
astronomical  literature  and  little  in  the  way  of  technical  guidance  (cf.  Young  1074, 
for  a  litany  of  photometry  pit  falls  relevant  to  this  program). 
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Wr  have  previously  found  (e.g.,  Lockwood  l‘)S-l )  that  intermediate-band  pho¬ 
toelectric  photometry  yields  nightly  mean  differential  magnitudes  witli  a  standard 
deviation  of  about  0.003  mag  (0.3%).  If  there  are  no  systematic  errors.  8 -1C  nights 
should  theoretically  yield  seasonal  mean  values  with  a  precision  on  the  order  of  0.001 
mag  (0.1%  ).  Obviously,  at  this  demanding  level  of  precision,  the  long-term  stability 
of  the  photometer  and  the  interference  filters  is  of  paramount  importance.  Ear¬ 
lier  long-term  monitoring  of  solar  system  objects  using  the  same  equipment  yielded 
results  that  appear  to  be  stable  to  about  0.003  mag  and  perhaps  better,  but  this 
prior  experience  offers  no  certa:n  proof  of  performance  at  the  even  higher  precision 
required  here  (Lockwood  ft  a!.  1980;  Lockwood  and  Thompson  1986u). 

Ultimately,  the  uncertainty  of  differential  magnitudes  propagates  from  the  mea¬ 
surement  errors  of  the  simplest  data  element,  the  photon-counting  statistics  of  indi¬ 
vidual  10-second  integrations,  plus  random  or  systematic  measurement  errors  due  to 
a  variety  of  causes.  As  a  fundamental,  observationally  testable  statistic  of  the  error 
budget,  we  estimate  how  precisely  the  differential  magnitudes  should  repeat  from 
cycle  to  cycle  within  the  nightly  sequence  y<b,b,y.  Known  or  suspected  components 
of  error  include  (1)  photon  counting  statistics  and  the  correction  for  dead  time, 
(2)  sky  background  measurements,  (3)  differential  extinction,  (4)  sky  transparency 
changes  during  the  cycle,  (5)  photometer  sensitivity  changes  during  the  cycle,  (6) 
scintillation,  and  (7)  effects  of  seeing  and  image  motion  within  the  star  diaphragm. 
In  the  following  sections,  we  estimate  as  best  we  can  the  relative  importance  of  these 
various  components. 

2.4.1  Photon-counting  errors 

The  throughput  of  the  telescope  +  photometer  system  changes  only  slowly 
with  time,  as  evidenced  by  the  constancy  of  stellar  and  Sr90  count  rates.  Photon- 
counting  errors  are  thus  a  simple  function  of  stellar  magnitude,  since  the  observations 
are  deliberately  made  in  an  identical  way  each  night.  The  stars  in  Table  1  have  a 
median  V  magnitude  of  6.1  mag;  the  faintest  is  8.2  mag.  Three-quarters  of  the  stars 
are  brighter  than  V  —  6.7  mag,  corresponding  to  an  observed  count  rate  in  the  y 
filter  >  25,000  sec'1,  and  greater  still  in  b  for  stars  hotter  than  B  -  V  ~  1.4  (K7V, 
K3II1).  Six  10-second  integrations  therefore  yield  an  internal  photon-counting  error 
<0.001  mag  for  V'  <  6.7  (i.e.,  for  75%  of  the  stars  in  Table  1).  In  §3.7.1  we  illustrate 
the  predicted  photon-counting  error  distribution  for  the  actual  set  of  differential  star 
pairs  in  the  observing  program. 

2. 1.2  Sky  background 

Under  dark  sky  conditions,  the  observed  count  rate  through  the  49  arcsec 
diameter  star  diaphragm  normally  used  was  ~250  sec-1 ,  yielding  a  negligible  (<1%) 
sky  contribution  to  the  stellar  signal  and  a  completely  negligible  ( <0.1%)  uncertainty 
in  the  resulting  sky-subtracted  count. 

An  unseen  12th  magnitude  background  star  will  contaminate  the  measurement 
of  a  7th  magnitude  star  at  the  1%  level,  whether  it  is  included  in  the  diaphragm 
with  the  star  or  is  accidently  encountered  during  the  separate  measurement  of  sky 
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brightness.  Fortunately,  the  density  of  12th  magnitude  stars  averaged  over  the 
whole  sky  is  fairly  low  (Roach  and  Gordon  1973):  the  probability  of  accidently 
including  a  background  star  in  a  49  arcsec  diaphragm  is  ~0.03  in  the  magnitude 
range  12.5  <  \  <  13.5  and  <0  05  for  all  stars  brighter  than  \  =  14.0,  a  negligible 
risk.  Nonetheless,  great  care  has  been  taken  to  avoid  including  faint  background 
stars  (see  §3.2). 

2.4.3  Extinction  and  airmass 

Scintillation  noise  and  atmospheric  extinction  effects  are  minimized  in  this  pro¬ 
gram  because  all  but  one  of  the  star  groups  transit  the  meridian  at  an  airmass  <1.3, 
Differential  airmasses  within  groups  are  usually  negligibly  small  (Figure  2).  If,  for 
example,  the  differential  airmass  range  is  as  large  as  0.03~  a  condition  occurring 
only  rarely— and  the  assumed  mean  extinction  coefficient  is  incorrect  by  0.03  mag 
airmass-1.  a  typical  uncertainty  (Lockwood  and  Thompson  198Gh),  the  resulting 
differential  magnitude  error  is  <0.001  mag.  Note,  however,  that  this  error  system¬ 
atically  affects  all  observations  of  a  particular  group  in  the  same  way  on  a  given 
night;  it  contributes  therefore  to  the  night-to-night  dispersion  and  can  in  principle 
masquerade  as  stellar  variability.  In  §3.7.2  we  test,  some  data  with  especially  “large” 
values  of  the  differential  airmass  for  possible  extinct  ion- related  systematic  effects. 

2.4.4  Temporal  atmospheric  transparency  changes 

The  observation  of  a  single  cycle  requires  at  most  10  minutes  for  a  quartet 
and  7  minutes  for  a  trio.  If  the  sky  transparency  changes  during  this  interval  (due 
to  vanishingly  thin  cirrus  clouds,  for  example),  random  errors  are  introduced  into 
the  differential  magnitudes.  Termination  of  observations  is  at  the  discretion  of  the 
observer  (see  §3.2);  but  generally,  problems  with  sky  transparency  become  noticeable 
only  when  they  reach  t lie  1-2%  level. 

Certainly,  on  truly  “photometric”  nights,  repeated  observations  of  bright  stars 
should  approach  the  scintillation  +  photon  limit.  In  some  earlier  measurements 
(cf.  Lockwood  1984),  we  looked  for  evidence  of  sky  transparency  fluctuations  over 
intervals  of  a  few  minutes  by  examining  the  relative  precision  of  the  differential 
magnitudes  of  the  various  pairs  of  stars  within  cycles.  Over  the  2-7  minute  time 
interval  in  which  the  various  pairs  of  stars  within  quartets  were  measured,  we  found 
no  evidence  for  systematic  effects — the  precision  attained  for  each  of  the  pairs  was 
the  same. 

2.4.5  Temporal  changes  in  photometer  sensitivity 

The  recorded  output  of  the  Sr90  standard  source  light  that  precedes  and  con¬ 
cludes  each  four-cycle  data  set  provides  confirmation  of  the  temporal  stability  of 
the  photometer  electronics.  Apart  from  a  slow  (<1%  hr-1)  drift  in  the  count  rate, 
evidently  related  to  ambient  temperature  changes  as  the  night  progresses,  the  source 
output  repeats  in  accord  with  photon-counting  statistics,  i.e.,  about  0.16%  rms  (10- 
second  integrations). 
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2. -1.6  Scintillation  error 


The  relevant  formulas  are  given  by  Young  (1974).  For  a  0.5-m  telescope,  the 
expected  standard  deviation  for  a  single  10-second  integration  is  0.0008  0.0018  mag 
in  the  aii  mass  range  of  interest,  depending  on  the  angle  between  the  line  of  sight  and 
the  upper-air  wind  direction.  Combining  six  10-second  integrations  thus  is  expected 
to  provide  an  essentially  negligible  0.0005  mag  contribution  to  the  error  budget. 

2.4.7  Errors  caused  by  seeing  and  image  motion 

The  observed  brightness  of  decentered  stars  can  vary  because  of  inhomogene¬ 
ities  in  the  sensitivity  across  the  focal  plane  aperture  of  the  photometer.  In  addition, 
decentered  or  moving  images  result  in  a  position-dependent  occultation  of  differ¬ 
ent  parts  of  the  extended  faint  wings  of  stellar  images  by  the  circular  diaphgram. 
Variable  seeing,  image  motion,  and  telescope  tracking  errors  combine  to  degrade 
measurement  precision  unpredictably. 

Drift  scans  in  perpendicular  directions  across  the  29-  and  19-arcsec  apertures 
indicate  that  the  photometer  response  is  uniform  within  0.3l/7  or  so  over  the  cen¬ 
tral  third  of  each  aperture  where  most  observations  are  made  (cf.  §3.2).  Random 
decentering  of  star  images  may  contribute  significantly  to  the  error  budget,  but  we 
can  only  guess  an  approximate  amount. 

2.4.8  A  hypothetical  error  budget 

A  hypothetical  error  budget  for  a  single  stellar  magnitude  determined  from  six 
10-second  integrations  includes  the  terms  in  Table  2.  In  this  example,  more  than 
half  the  variance  is  conceivably  contributed  by  the  poorly  estimated  component 
related  to  image  motion.  When  two  individual  stellar  measurements  are  combined 
to  create  a  differential  magnitude,  the  variance  is  doubled,  leading  to  an  expected 
cycle-to-cycle  rms  dispersion  of  0.002-0.005  mag.  In  §3.7.1,  we  compare  the  error 
budget  with  real  data,  and  in  §3.7.2  we  will  show  that  the  value  0.002  corresponding 
to  the  low  estimate  above  is  achieved  for  the  most  stable  stars. 

3.  THE  DATA 

3.1  Observations  obtained 

The  stars  listed  in  Table  1  are  organized  into  30  groups  of  trios  and  quartets. 
Trios  were  promoted  to  quartets  if,  during  the  course  of  the  program,  more  than 
one  star  was  found  to  vary.  Nearly  2000  nightly  observations,  about  00  nights  per 
group,  were  obtained  on  350+  nights  between  March  1984  and  December  1987  — 
three  full  seasons  for  all  groups  and  four,  or  part  of  a  fourth,  for  nearly  all  groups. 
All  observations  were  made  and  reduced  to  differential  magnitudes  by  B.  A.  Skiff. 
A  t>  pical  original  data  record  is  given  in  the  Appendix. 

3.2  Observational  procedures 

Observations  were  about  evenly  split  between  diaphragm  3  (29  arcsec)  and 
diaphragm  4  (49  arcsec).  One  or  the  other  was  selected  based  on  the  apparent 
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TABLE  2.  Error  Budget 


Low 

Estimation 

(mag) 

High 

Estimation 

(mag) 

Photon  counting  (V  =  6.7) 

0.001 

0.001 

Sky  subtraction 

0.0005 

0.001 

Differential  extinction 

<0.0005 

0.001 

Scintillation 

0.0003 

0.0007 

Sensitivity  drift 

.  -  . 

0.0005 

Centering/ 1  racking 

0.001 

0.003 

Root  mean  square 

0.0016 

0.0036 

image  size  of  bright  ( l ’  <  2)  stars:  if  the  seeing  disk  was  less  than  6  to  8  arcsec 
(corresponding  roughly  to  FVVHM  and  estimated  using  the  known  sizes  of  the  di¬ 
aphragms),  then  the  smaller  aperture  was  used;  worse  seeing  dictated  the  use  of  the 
larger  aperture.  Observations  were  not  made  (or  halted  at  some  convenient  point) 
if  the  images  were  larger  than  about  25  arcsec,  a  condition  that  occurred  a  few 
nights  each  year  associated  with  passage  of  winter  storms.  The  estimated  seeing  is 
recorded  in  the  comments  line  cf  each  file  header. 

Observations  were  begun  during  late  twilight  when  the  sky  counts  declined 
to  about  double  the  expected  full-dark  level  for  the  appropriate  diaphragm  size. 
Measurements  were  not  made  closer  than  20  to  30  degrees  from  the  Moon  or  when 
the  sky  counts  exceeded  1000  sec-1  in  the  b  filter  due  to  moonlight.  At  high  galactic 
latitudes,  the  sky  readings  were  made  by  simply  flipping  the  declination  slow-motion 
crank  either  north  or  south  about  a  full  turn;  tests  showed  this  corresponds  to  a 
typical  offset  "throw”  of  -4  or  5  arcminutes.  Except  during  the  first  few  months  of 
t he  project,  the  first  two  cycles  have  sky  readings  to  the  south  and  the  last  two  to 
the  north.  At  low  galactic  latitudes,  sky  readings  were  taken  at  locations  selected 
during  dark  time  and  used  consistently;  for  several  groups  this  meant  taking  all  the 
sky  readings  in  one  direction.  All  the  low-latitude  groups  could  be  measured  using 
diaphragm  4. 

Stars  were  carefully  centered  in  the  diaphragm,  usually  within  the  central  one- 
third  of  the  diameter.  The  telescope  drive  often  exhibited  a  slow  drift  caused  by  a 
small  periodic  error  in  the  worm,  and  stars  were  usually  offset  to  one  side  or  the 
other  of  the  central  one-third  of  the  field  to  let  the  star  drift  across  dead  center.  For 
the  numerous  double  stars  where  both  components  were  visible,  pairs  were  set  on  the 
weighted  photocenter  (determined  roughly  by  eye),  which  generally  was  very  close 
to  the  primary  star.  The  fields  were  usually  observed  within  one  hour  of  meridian 
passage,  and  rarely  up  to  2.5  hours  hour  angle. 
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Contamination  of  both  star  and  sky  founts  by  faint  field  stars  was  anticipated 
both  in  planning  and  in  observing.  Bright  stars  were  thoroughly  examined  by  var¬ 
ious  double  star  observers  about  100  years  ago  with  large  telescopes.  Thus,  any 
significant  close  companions  are  recorded  in  double  star  catalogs  (and  noted  in  Ta¬ 
ble  1).  Many  more  of  our  program  stars  are  also  in  double  star  catalogs  by  virtue 
of  their  large  proper  motions  micrometric  measurements  of  faint  background  stars 
being  used  to  ascertain  this  independent  of  meridian  circle  measurements.  These 
ended  up  being  listed  as  double  stars  because  not  a  few  of  the  faint  companions 
turned  out  to  have  common  motion.  At  any  rate,  the  point  is  that  we  know  about 
reasonably  bright  (V  <  11)  stars  near  the  progrant/comparison  stars  that  would  be 
included  in  the  diaphragm. 

The  declination  offset  throws  were  not  done  completely  blindly:  the  “WP1IOT” 
computer  program  that  records  data  warns  the  observer  of  overbright  readings  on 
the  sky  (>  ‘2x  normal).  The  readings  that  were  noted  as  even  slightly  above  normal 
(indicated  by  the  one-second  integrations  on  a  frequency  counter)  were  repeated  at 
a  slightly  different  location  to  try  to  get  the  offending  star  out  of  the  diaphragm  as 
the  integrations  continued.  'Lgarding  the  discussion  ip.  §2.4  ‘2.  stars  V  <  14  were 
hardly  ever  included  in  sky  readings,  not  so  much  because  the  odds  of  inclusion  are 
so  small,  but  because  the  observing  technique  specifically  attended  to  this  point. 
Certain  stars  have  a  faint  ••companion”  consistently  falling  in  the  normal  “sky” 
position,  and  these  have  all  their  sky  readings  taken  in  the  opposite  direction. 

Finally,  it  is  also  worth  noting  that  the  sky  background  at  low  galactic  latitudes 
is  noticeably  brighter  due  to  the  presence  of  stars  below  the  limit  of  visibility  of  the 
telescope  ( V  >  15  on  dark  nights). 

3.3  Data  rejection  criteria 

Integrations  were  monitored  regularly  for  anomalous  counts  both  as  they  were 
taken  and  using  a.  summary  facility  in  the  photometry  program.  Generally,  the 
batches  of  repeated  10-second  “star”  integrations  would  show  a  total  range  of  1%  or 
a  bit  more;  larger  ranges  were  indicative  of  clouds,  miscentered  stars,  or  poor  seeing 
(i.o.,  too-small  a  diaphragm).  A  set  of  integrations  would  be  repeated  only  due  to 
observer  or  mechanical  error,  viz.  the  wrong  star,  miscentering  due  to  drive  motor 
drift  or  wind.  The  rms  noise  on  a  set  of  six  integrations  was  typically  0.3-0. 5%;  on 
the  most  stable  nights  (generally  associated  with  the  best  seeing),  the  rms  would 
consistently  reach  0.1  0.3%;  and  on  poor  nights,  increase  to  0.5  0.8% — these  last 
often  preceding  the  arrival  of  cirrus,  easily  discernable  photometrically  at  the  1-2% 
level.  The  repeated  cyclic  measurements  provided  further  checks  on  photometric 
stability.  If  a  star  did  not  repeat  to  within  about  half  a  percent  between  cycles 
(after  compensating  mentally  for  airmass  changes),  suspicions  were  raised;  usually, 
continued  observations  would  show  persistent  deviations  (either  cycle-to-cycle  or 
intra  cycle),  thus  ending  the  night. 


.'{.  I  Reductions 


Utilizing  data  reduction  procedures  that  have  been  highly  standardized  for 
many  years,  the  raw  observations  were  corrected  for  pulse  pair  “dead  time.”  the  sky 
background  was  subtracted,  and  differential  magnitudes  were  computed  using  sea¬ 
sonally  adjusted  mean  values  of  the  extinction  coefficient  ■  (Lockwood  and  Thompson 
19866).  The  principal  raw  data  archive  is  maintained  on  the  original  8-inch  flexible 
disks. 


Each  reduced  observation  yields  a  printed  output  and  one  data  “record”  per 
cycle  that  is  appended  to  the  master  file  for  each  group.  Every  cycle  produces  a 

set  of  three  differential  magnitudes  for  a  trio  (star  1-star  2,  1-3,  2-3),  six  for  a 

quartet  (star  1-star  2.  1-3,  1-4,  2-3,  2—4,  3-4).  Examples  of  the  printed  nightly 
output  file  are  shown  in  the  Appendix. 

3.5  Seasonal  mean  differential  magnitudes 

Table  3  contains  the  seasonal  mean  Julian  Date,  the  mean  differential  b  and 
y  magnitudes,  and  their  standard  deviations  for  each  group.  Nomenclature  of  the 
pairwise  magnitude  combinations  corresponds  to  the  star  numbers  assigned  within 
each  group  in  Table  1  The  number  of  observations,  rz,  is  the  number  of  cycles 

obtained  in  each  filter:  hence  the  number  of  nights  is  n/2. 

3.6  Light  Curves 

Figures  4.1-4.30  show  the  light  curves  for  each  group.  Because  the  light  curves 
in  b  and  y  are  essentially  identical,  the  data  have  been  condensed  for  presentation 
by  averaging  the  y  magnitudes  of  cycle  1  with  the  b  magnitudes  of  cycle  2,  and  the 
b  magnitudes  of  cycle  3  with  the  y  magnitudes  rf  cycle  4.  Individual  data  points 
for  the  averaged  ( b  -f  y)/ 2  cycles  are  shown  on  the  left-hand  panels,  using  the  same 
vortical  scale  (0.05,  0.1,  or  0.2  mag  full  range  as  required)  for  all  the  panels  of  a 
given  group. 

Seasonal  mean  magnitudes  are  summarized  on  the  right-hand  panels,  where 
the  mean  values  are  connected  by  solid  lines.  The  error  bars  are  95%  confidence 
intervals  about  the  means.  Median  differential  magnitudes  are  shown  as  open  cir¬ 
cles:  by  comparing  the  relative  locations  of  means  and  medians,  some  feeling  for  the 
asymmetry  of  the  data  can  be  obtained,  particularly  useful  in  assessing  the  impor¬ 
tance  of  outliers  that  may  be  observed  on  the  left-hand  panels.  Note  that  the  95% 
confidence  intervals  are  sensitive  to  not  only  the  rms  dispersion  of  the  data,  but  also 
to  the  number  of  data  points:  large  confidence  intervals  often  indicate  sparseness  of 
data  rather  than  variability. 

3.7  Observed  precision  levels 

Before  attempting  to  determine  which  stars  may  bo  varying,  it  is  important  to 
establish  a  clear  picture  of  the  magnitude  of  various  error  sources;  i.e.,  a  statistical 
profile  of  the  data  as  a  whole.  In  this  section  we  partition  the  data  in  different  ways 
to  examine  (1)  cycle-to-rycle  magnitude  variations,  (2)  errors  related  to  differentia] 
extinction,  and  (3)  errors  associated  with  bad  nights  or  a  particular  star  group. 
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TABLE  3.  Summary  of  Annual  Mean  Differential  Magnitudes 


o 

o 

o 


o  o  o  o 
oooo 
oooo 


O  -  X!  05 
cD  r-  cd  uo 
04  04  04  04 
«  M  cc 


O  0  —  1-04  go  —  O  CD  iO  —  X 

o  oooo  oooo  oooo 

o  oooo  oooo  oooo 

o  oooo  oooo  oooo 


x  —  x  x  x  04  cd  i- 

O  04  —  O  I  —  l—  cd  cd 

x  x  x  x  to  iO  tc  to 

oooo  oooo 
i  I  i  I 


-r  —  x  co  cd  i.o  re  cd 

CM  04  —  —  —  CD  CD  CD 

o  q  o  o  i-  r-  o  io 

x  x  x  x  oc  x  x  x 

O  o  o  o  —  — '  —  — 

1111  till 


oooo 

oooo 


lO  rr  lO  X) 
O  X  —  04 
OC  (—  CD 

04  04  04  04 

1  I  i  I 


OOOO 

oooo 

oooo 


r—  cc  — i  r— 

lO  lD  X  04 
04  04  04  04 
C3  C3  03  os 

OOOO 

I  I  1  I 


O  O  04  o 

oooo 

oooo 


O  x  —  o 

CD  to  X  rr 

x  x  x  x 
oc  oc  OO  CO 

oooo 

1111 


x  04  —  -r 
x  x  x  x 
oooo 
oooo 


r-  o  o  x 

CO  04  CO  04 
OOOO 
OOOO 


OOOO  OOOO 


«S  Q  O 

—  2 


r—  ^  t-  o 

CO  I  —  CD  CD 
X  C3  X  03 

oooo 
oooo 
i  I  (  I 


X  —  O  — 
CO  CO  CO  CO 

oooo 

oooo 

oooo 


r—  to  id  o 
o  O’  CO  — 
O  04  — i  04 
—  O’  "T  O' 

oooo 

till 


OC-  CO  X  X 
04  lO  CO  04 


O  N  50  tO 

co  ~r  — •  x 
O  CO  h  O 
CD  CD  CD  T— 


ic  d?  ift 
iO  CD  X  iO 


OOOO 

I  I  I  I 


lO  I'-  —  03 
CO  04  CO  04 
OOOO 
OOOO 


CO  LO  f—  CO 
CD  X  CC  CO 
t-  x  r-  x 
rr  rr  rr  •'T 
oooo 
I  I  I  I 


X  co  X  X 
O 1  ‘O  X  04 


cd  *-r  *-r  co 
©  t—  co  »c 
CO  •O’  —  X 
0X1-0 
CD  CD  CD  C  — 


OOOO 

OOOO 


O  O  CD  CO 
05  O  O  >0 
O  —  OO 
lO  tO  tO  tO 


X  o  X  o 
0  —  00 
oooo 
oooo 
oooo 


*D  Q  to  o 

03  ©  X  CO 


OOOO 

OOOO 


CD  CM  O  — 
—  —  —  CM 
OOOO 
OOOO 


OOOO  oooo 


X  X>  O  X  OO  to  'T  05 

—  —  —  —  o  —  —  — 
oooo  oooo 
oooo  oooo 

OOOO  OOCDO 


tilt  lilt 


cd  X  o  r- 
0  0  —  0 
OOOO 
oooo 


CD  t-  —  o 

04  f  1  X  05 
O  OX  X 

oooo 

6666 

till 


O  O  X  X 
X  X  04  CM 


CO  "T  C  4  O 
—  CD  —  X 
0X1-0 
CD  CD  CD  I— 


—  O  X  lO 

—  —  o  o 
oooo 
oooo 


O  CD  04  — 

t-  cd  -r  --r 
t-  t  ’? 
OOOO 


O  O  X  X 
X  X  04  04 


x  x  — 

CD  *r  04  o 
—  CD  —  X 
O  X  t—  o 
co  CD  CD  I— 


O  O  rr  o 

04  CM  —  — 
iO  tO  tO  lO 

t-  r-  l—  i- 
6  6  6  6 
till 


CM  X  X 

—  —  o  — I 
oooo 
oooo 


CD  lO  O  t— 
CM  X  X 
CM  CM  —  O 


rr  x  o  rr 

04  —  04  — 


O  —  —  X 
O  X  X  o 
o  x  r—  o 

CD  CD  CD  t- 


X  O  t-  t- 
X  t—  rr  rr 

oooo 


CO  CD  —  X 
t— *  CD  •— h  »— i 

oooo 

oooo 


CD  X  to  X 
X  X  X  X 
x  x  x  x 
oooo 
till 


oooo 

oooo 


X  LO  —  O 
X  CM  CO  o 
03  03  X  t— 
C4  CM  CM  CM 


rr  x  o  -r 

04  —  CM  — 


©  ©  03  X  XCDXX 

03  —  —  X  iC  03  03  C3 

OXXC3  —  x  X  o 

o  x  r-  o  o  x  r-  o 

co  CO  co  t~—  CD  CD  CD  t— 


03  03  X  — 

ci  ri  o  05 
rr  rr  rr  x 
OOOO 


lit!  I  1  I  I 


OOOO 

OOOO 


oooo  oooo 


CD  X  LO  o 
CM  CM  «— *  CM 
X  X  X  X 

oooo 
I  I  I  I 


03  —  X  — 
CD  t-  X  iO 
OOOO 
CM  CM  CM  CM 

oooo 
I  I  I  I 


»  O  113  H 


X  CD  X  X 
Ut>  03  03  03 
—  X  X  03 
O  X  t-  O 
CD  CD  CD  t— 
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TABLE  3.  Summary  of  Annual  Mean  Differential  Magnitudes  (continued) 
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5841.3  18  -  0.9408  0.0023  -  0.8519  0.0029  -  0.5031  0.0047  0  0889  0.0018  0.4377  0.0040  0.3488  0.0030 

0223.4  20  -0.9428  0.0027  -0.8493  0.0033  -0.5007  0.0030  0.0935  0.0024  0.1301  0.0021  0.3427  0.0028 

6011.5  20  -0.9501  0.0023  -0.8729  0.0036  —0.5123  0.0034  0.0772  0.0031  0.4378  0.0027  0.3000  0.0038 

6965.4  20  -0.9400  0.0017  -0.8675  0.0031  -0.5043  0  0025  0.0725  0.0035  0.4357  0.0024  0.3031  0.0045 
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1  ABI.K  3.  Summary  of  Annual  M*>an  Differential  Magnitudes  (continued) 
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Figure  4.1.  Light  curves  of  the  HD1835  group,  compressed  scale,  b  and  y  averaged,  (left)  individual  data  po'  .ha 
y  averaged,  (right)  seasonal  mean  magnitudes  (connected  by  solid  lines),  and  medians  (solid  dots).  Krror  bars  denote  9' 
confidence  intervals.  Star  1  =  111)1835  (G2V)  is  strongly  variable  and  has  an  observed  photometric  rotation  period  of  7.7  d? 
(Table  6). 


Figure  4.2.  Light  curves  of  the  HD10476  group.  Star  1  =  HD10476  (K1V)  is  variable  from  year  to  year,  but  all  three  stars  are 
unusually  stable  each  season.  The  range  of  star  1-star  3  is  only  0.0009  mag  in  four  seasons. 


ZZ I 


Julian  Date  Julian  Date 

Figure  4.3.  Light  curves  of  the  1ID13421  group.  Star  2  =  HD13611  (C(iII  III),  an  MK  standard  star,  is  strongly  variable  from 
year  to  year. 


ire  4.5.  Light  curves  of  the  HD23140  group.  Star  1  =  HD23140  (K2)  varied  in  two  seasons,  but  none  of  the  stars  showed 
term  variability. 


-1.4 

1.42 


curves  of  the  HD35296/HD39587  group.  Star  1  =  HD35296  (F8V)  and  star  3 
?d  strongly  from  year  to  year  as  well.  We  find  photometric  rotation  periods 
agreement  with  the  II K  period.  One  or  both  of  the  comparison  stars  appe 
ear  to  year.  Note  the  variation  of  star  2  just  at  tiie  end  of  the  second  season 


Fig  -re  4.8.  Light  curves  of  the  HD76572/HD75332  group.  All  three  of  the  original  trio  pairs  are  very  stable,  with  the  four-season 
range  of  star  2-star  3,  0.0007  mag,  the  smallest  observed.  Star  2  =  HD73596  (F5II1)  may  vary  slightly  on  short  time  scales. 
Star  4  =  HD75332  (F7Vn)  is  very  active. 


Figure  4.9.  Light  curves  of  the  HD81809  group.  Comparison  star  3  =  HD82074  (G6IV)  varied  strongly  from  year  to  year.  Star 
=  HD81809  (G2V)  may  fluctuate  slightly. 


Figure  4.11.  Light  curves  of  the  HD  103095  group.  Star  3  =  HD101606  (F4V),  a  why  standard,  may  be  slightly  variable,  long 
term,  but  this  group  is  essentially  quiescent. 


ie  HD114710  group.  Star  3  =  HD112980  (G9III)  is  consistently  and  strongly  variable. 


Figure  4.13.  Light  curves  of  the  HD115383/IID117176  group.  Star  1  =  HD1 15383  (GOVs),  a  uvby  and  MK  standard  star,  varies 
strongly  from  year  to  year. 


Figure  4.14.  Light  curves  of  the  HD115404  group.  Star  1  =  HD115404  (K1V  +  M1V)  varies  strongly.  The  derived  photometric 
rotation  period  of  11.0  days  is  probably  spurious:  it  is  not  in  accord  with  the  well-determined  HK  period  of  18.8  days.  Star  ‘2-star 
3,  though  quiescent,  vary  from  year  to  year. 


Figure  4.15.  Light  curves  of  the  HD120136  group.  Star  1  =  HD  120 136  (F6IV)  seems  to  vary  strongly  from  year  to  year  but 
neither  of  the  comparison  stars  is  particularly  stable. 


Figure  4.16.  Light  curves  of  the  HD124570  group.  Star  3  =  HD126246  (F8V  4-  G1V)  varies  strongly  with  a  range  of  nearly 
0.02  mag. 


■  4.17.  Light  curves  of  the  HD129333  group,  compressed  scale.  Star  1  =  HD129333  (GO)  has  the  largest  log  v; 
program  and  the  greatest  range  of  variability  (0.05  mag).  Its  photometric  rotation  period  in  two  seasons  was  2.7  and 


ght  curves  of  the  HD143761  group.  Star  1  =  HD143761  (GO+Va),  a  uvby ,  MK,  and  HK  standard  is  exceptionally 
HD140716  (gG9)  may  fluctuate  slightly. 


vox- 


Figure  4.20.  Light  curves  of  the  HD149661/HD152391  group,  compressed  scale.  Star  1  =  HD149661  (K2V)  and  star  3  = 
HD152391  (G7V)  both  vary  strongly  with  ranges  of  0.01  and  0.02  mag,  respectively.  One  or  both  comparison  stars  is  slightly 
active  also,  but  the  long-term  range  is  small,  about  0.003. 


Figure  4.21.  Light  curves  of  the  HD157856/HD158614  group.  Though  consistently  quiescent,  star  1  =  HD157856  (F3V)  may 
show  long-term  variation.  Star  3  =  HD158614  (G8IV)  varies  strongly  from  year  to  year.  Both  comparison  stars  are  stable. 


Light  curves  of  the  HD160346  group.  This  group  lacks  a  quiescent  pair  of  comparison  stars.  Star  2  =  HD160385 
variable  every  year  but  is  stable  long  term.  Star  3  =  HD160823  (GO)  varies  strongly  from  year  to  year. 
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Figure  4.23.  Light  curves  of  the  HD161239  group.  Both  star  2  =  HD162211  (K2III)  and  star  3  =  HD  162076  (G51V)  fluctuate 
slightly,  and  star  4  =  1ID 160935  (F5)  seems  to  vary  strongly  from  year  to  year.  The  long-term  light  curves  are  confused. 


Figure  4.25.  Light  curves  of  the  HD182572  group.  Star  3  =  IID185018  (GOIb-11)  may  vary  slightly  from  year  to  yea: 


Light  curves  of  the  HD190007  group.  Star  1  =  JID190007  (K4)  varies  strongly  on  all  lime  scales. 


:.28.  Light  curves  of  the  HD201091/HD201092  group.  In  this  quiescent  group,  only  star  1  =  11D201091  (K5V),  a  uvby 
tandard.  is  consistently  variable  each  year.  Its  companion,  star  3  =  HD201092  may  fluctuate  slightly.  Star  1  and  star  1 
from  season  to  season,  but  the  light  curves  are  confused. 


of  the  HD215704  group,  compressed  scale.  Star  2  =  HD215427  (K5)  varies  s 


HD216385  group.  Star  1  _  IID216385  (F7IV)  varies  from  season  to 


3.7.1  Cyc’^-to-cycle  differences 

Comparing  the  differential  magnitudes  of  the  same  pair  of  stars  observed  in 
two  separate  cycles  on  the  same  night  provides  a  useful  measure  of  the  precision  ac¬ 
tually  attained  on  ~l/2-hour  time  scales,  where  the  inconstancy  of  the  atmospheric 
transparency  and  photometer  response  may  be  possible  sources  of  noise.  On  this 
time  scale,  all  the  stars  are  assumed  to  be  nonvariable. 

Because  the  cycle  order  is  always  y,b,b,y,  the  comparison  of  cycles  2  and  3 
(separated  by  less  than  10  minutes)  and  cycles  1  and  -1  (separated  by  about  20 
minutes)  provides  a  diagnostic  of  the  stability  of  measurement  over  these  intervals. 
We  computed  the  absolute  values  of  the  differences  [cycle  1 -cycle  4]  and  [cycle 
2-cycle  3]  to  investigate  the  repeatability  of  measurements  independent  of  night- 
to-night  intrinsic  stellar  variability  and  measurement  error.  Because  this  quantity 
is  also  independent  of  the  actual  differential  magnitudes  of  the  stars,  it  can  be 
summed  over  all  groups  and  all  nights  to  provide  global  statistics  pertaining  to 
the  entire  program,  not  just  small  data  subsets.  Note  that  the  absolute  cycle-to- 
cycle  differences  are  larger  by  a  factor  \/2  than  the  rms  cycle-to-cycle  dispersion  of 
differential  magnitudes. 

The  top  panel  of  Figure  5  shows  the  observed  distribution  of  the  absolute 
values  of  cycle-to-cycle  differences  of  differential  magnitudes  for  a  representative 
subset  of  all  the  observations.  The  observed  error  distributions  are  similar  in  shape, 
with  medians  of  0.0024  and  0.0022  and  upper  quartiles  0.0044  and  0.0039  mag  for 
the  absolute  values  [cycle  1 -cycle  4]  ( y )  and  [cycle  2-cycle  3]  (6),  respectively. 
Higher  count  rates  in  b  for  all  stars  earlier  than  K3  probably  accounts  for  the  slight 
differences  between  the  two  histograms  rather  than  the  temporal  proximity  of  the  b 
cycles. 

In  the  bottom  panel  of  Figure  5,  the  predicted  cycle-to-cycle  photon-counting 
error  distribution  is  shown,  based  on  the  actual  B  and  V  magnitudes— and  the  cor¬ 
responding  known  count  rate — of  the  66  pairs  of  stars  comprising  the  22  original  trio 
groups.  (A  distribution  computed  for  quartet  pairs  should  be  essentially  identical 
but  was  not  computed.)  Note  that  the  variances  contributing  to  this  distribution  are 
a  factor  of  four  larger  than  the  photon  count  variance  for  a  single  star:  computation 
of  a  differential  magnitude  doubles  the  variance,  and  subtracting  one  differential 
cycle  from  another  doubles  it  again. 

Comparing  the  two  panels  of  Figure  5.  we  observe,  with  no  surprise,  that  real 
observations  ( top  panel)  evidently  include  contributions  to  the  total  variance  that 
greatly  exceed  the  photon  counting  error  ( bottom  panel).  Some  guidance  to  the 
relative  importance  of  various  sources  of  error  is  suggested  in  Table  2;  but  as  the 
importance  of  the  largest  suspected  variance  source  (i.e.,  centering)  is  basically  just 
a  guess,  we  must  conclude  tentatively  that  it  is  the  small  atmospheric  transmission 
fluctuations  on  less-than-perfect  photometric  nights  that  contribute  most  to  the  tail 
of  the  histograms  of  Figure  5 a. 


Figure  5.  (top).  Observed  distribution  of  the  absolute  values  of  cycle-to-cycle 
differences  obtained  within  the  y.b,b,y  sequence  each  night.  A  representative  subset 
of  the  whole  data  set  has  been  used.  ( solid  line )  -  y  data.  ( dotted  line)  -  b  data. 
(bottom).  Predicted  cycle-to-cycle  error  distribution  due  only  to  photon  statistics, 
based  on  the  V  and  B  magnitudes  of  the  pairwise  combinations  of  stars  in  the  30 
groups.  ( solid  line)  -  V.  ( dotted  line)  -  B. 

3.7.2  Errors  arising  from  differential  extinction 

Figure  6  shows  the  actual  distribution  of  differential  airmasses  for  the  whole 
set  of  observations.  The  upper  quartile  of  the  distribution  has  the  value  0.02  air- 
mass.  Assuming  an  uncertainty  in  the  assumed  seasonally  adjusted  mean  extinction 
coefficient  ~0.03  mag  airmass-1  (Lockwood  and  Thompson  19866),  the  resulting 
differential  magnitude  error  is  ~0.0006  mag,  a  negligible  addition  to  the  error  bud¬ 
get  contributing  insignificantly  to  the  night-to-night  dispersion  of  nonvariable  stars. 

For  those  few  groups  having  “large”  (>0.02)  differential  airmasses,  averaging 
the  nightly  cycles  will  not  reduce  the  systematic  errors  arising  from  incorrect  as¬ 
sumed  values  of  differential  extinction.  Figure  7  shows  the  relation  between  mean 
intraseason  dispersion  and  differential  airmass  for  a  small  group  of  stars  presumed 
to  be  intrinsically  nonvariable  on  the  basis  of  their  low  chromospheric  activity  (log 
f?'HK)  values.  There  is  no  tendency  for  the  dispersion  to  increase  at  larger  values  of 
differential  airmass,  so  we  conclude  that  the  use  of  mean  extinction  coefficients  is 
not  a  significant  source  of  systematic  error. 
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Mean  rms  dispersion 


Figure  0.  Distribution  of  differential  airmasses  for  ail  of  the  observations.  The 
third  quartile  of  the  distribution  occurs  at  0.02  airmass. 


Figure  7.  The  relation  between  the  seasonal  rms  dispersion  of  differential  magni¬ 
tudes  and  differential  airmass  for  stars  unlikely  to  be  intrinsically  variable. 
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3.8  Possible  instrumental  effects 

3.8.1  Color  effects 

Over  the  length  of  an  observing  season  the  conditions  under  which  observations 
of  a  given  group  are  made  can  change  significantly.  For  example,  between  January 
and  June  the  mean  nighttime  low  temperature  rises  from  —  10°C  to  +5°C  (NOAA 
Local  Climatological  Data  for  Flagstaff,  Arizona).  To  test  for  seasonal  temperature- 
related  effects  or  temporal  effects  of  unknown  origin,  we  plotted  on  Figure  8  the 
cycle-to-cycle  rms  dispersion  ( b  and  y  averaged)  as  a  function  of  the  difference  in 
b  -  y  color  of  the  two  stars  of  a  pair  for  pairs  that  showed  no  evidence  of  intrinsic 
variability  in  any  season.  Temporal  (perhaps  temperature-related)  changes  in  the 
transmission  of  either  filter  should  cause  the  star  pairs  with  large  differences  mb  —  y 
color  to  exhibit  a  larger  dispersion.  No  effect  is  evident. 

3.8.2  Brightness  and  brightness  difference  effects 

For  the  brighter  stars,  dead-time  corrections  are  not  insignificant,  amounting  to 
more  than  0.5%  at  V  =  5.0.  Temporal — perhaps  temperature-related — changes  in 
the  amplifier  output  pulse  width  could  cause  systematic  effects.  Changes  in  linearity 
of  the  pulse-counting  electronics  will  be  manifested  as  an  increase  in  the  dispersion 
as  a  function  of  the  magnitude  difference  of  the  two  stars  of  a  pair  and  as  a  function 
of  the  magnitude  of  the  brighter  member  of  a  pair.  Using  the  same  data  as  in  §3.8.1 
above,  we  plotted  these  relationships  in  Figures  9  and  10.  No  effects  are  evident. 

1.  INTRINSIC  STELLAR  VARIABILITY  ON  INTRASEASON  TIME  SCALES 

Observations  in  just  one  season  provide  information  on  the  variability  of  stars 
on  time  scales  of  days  to  months  and  can  reveal  the  modulation  of  disk-integrated 
brightness  caused  by  the  flux  deficits  as  spotted  regions  rotate  into  view.  On  longer 
time  scales,  the  evolution  of  spot  groups  may  in  principle  be  inferred  by  changes  in 
the  amplitude,  phase,  and  period  of  rotational  light  curves.  Having  developed  the 
tools  needed  to  distinguish  low  level  intrinsic  stellar  variations  from  observational 
artifacts,  we  proceed  to  characterize  the  intrinsic  variations. 

When  a  group  contains  only  one  varying  star,  as  is  ideally  but  not  always  the 
case,  its  identification  is  uriarnbigous:  the  two  (three  for  a  quartet)  differential  light 
curves  for  star  pairs — each  including  the  variable  as  one  of  its  components — are 
correlated.  Variability  is  thus  quantitatively  confirmed  by  the  attained  significance 
level  of  the  correlation  coefficients  and  is  supported  usually  by  an  obvious  pattern 
in  the  rms  dispersions  of  the  several  differential  magnitudes. 

•1.1  Intraseason  standard  deviations 

The  rms  dispersion  of  the  seasonal  mean  differential  magnitudes,  and  sy, 
based  usually  on  two  observations  per  night  in  each  filter,  can  signal  the  presence  of 
intrinsic  stellar  variability.  For  a  given  star  pair,  .$(,  ~  so  to  reduce  tabulation  in 
this  report,  we  shall  use  the  standard  deviations  of  the  combined  b  and  y  magnitudes, 
,s t,y,  formed  by  averaging  the  y  magnitudes  of  cycle  1  with  the  b  magnitudes  of 
cycle  2,  and  the  b  magnitudes  of  cycle  3  with  the  y  magnitudes  of  cycle  4,  thus 
producing  two  nightly  data  points  for  each  group  observed. 
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Figure  8.  Average  intraseason  standard  deviations  (b  and  y  averaged)  as  a  function 
of  the  difference  in  (6  -  y)  color  for  pairs  of  stars  showing  no  evidence  of  intrinsic 
variability. 


Figure  9.  Average  intraseason  standard  deviations  ( b  and  y  averaged)  as  a  function 
of  the  difference  in  V  magnitude  for  pairs  of  stars  showing  no  evidence  of  intrinsic 
variability. 


Figure  10.  Average  intraseason  standard  deviations  ( b  and  y  averaged)  as  a  func¬ 
tion  of  the  V  magnitude  of  the  brighter  member  of  pairs  of  stars  showing  no  evidence 
of  intrinsic  variability. 
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(rms  disp.  b  mags)/(rms  disp.  y  mags) 


Table  4  gives  sby  values  for  each  observing  season,  the  Julian  Date  range,  and 
n,  the  number  of  averaged  [(6  -f  y)/2]  data  points.  Except  for  the  rare  occasions 
when  only  one  b  and  one  y  cycle  were  observed,  the  number  of  nights  each  season 
is  therefore  nf 2.  In  each  group,  star  1  provides  the  group  name  and  is  the  principal 
program  star;  the  second  program  star  in  the  group  label  for  certain  quartets  is 
star  3  unless  otherwise  noted. 

Following  each  sby  value  in  Table  4  is  an  indication  of  whether  that  particular 
star  pair  shows  variability:  V  -  definitely  variable  (99%  confidence),  v  -  possibly 
variable  (95%  confidence),  or  C  -  constant,  according  to  the  significance  of  corre¬ 
lation  coefficients  (§4.3).  These  designations  are  omitted  when  the  identity  of  the 
variable  star(s)  is  ambiguous. 

4.2  The  relation  between  sb  and  sy 

Although  sb  ~  sy,  the  ratio  sb/sy  ~  1.1  for  variables  and  ~  1.0  for  constant 
stars,  reflecting,  at  large  values  of  sb  and  sy,  mainly  an  aspect  of  intrinsic  variability 
and,  at  small  ( <0.004  mag)  values  of  sb  and  sy,  mainly  observational  noise.  Figure  11 
shows  sb/sy  plotted  as  a  function  of  sb.  For  stars  whose  consistent  fluctuations 
every  year  verified  their  intrinsic  variability,  we  find  sb  ~  1.15sy  and  for  consistently 
nonvariable  stars,  sb  ~  sy. 


0  .005  .01  .015  .OS 

rms  dispersion  of  b  magnitudes 

Figure  11.  The  ratio  of  the  intraseason  standard  deviation  in  b  to  the  standard 
deviation  in  y,  sb/sy  as  a  function  of  sb.  (filled  circles)  ■  consistently  constant  stars. 
(open  circles)  -  consistently  variable  stars.  The  ratio  is  typically  ~1.1  for  variable 
stars. 
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TABLED.  Intraseason  Dispersion  and  Interseason  Ranges  in  (b+y)/2 
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TABLE  4.  lntraseason  Disporsion  ami  Interseason  Hanses  in  (b+y)/2  (conf  intKii) 
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f  Star  4  =  HD75332 
;  Star  2  =  HD82635 


TABLE  4.  Intraseason  Dispersion  and  Interseason  Ranges  in  ( b+y)/ 2  (continuer!) 


HD143761  22  5782.9-5866.7  0.0021  0.0033  0.0026 

18  6145.0-6222.8  0.0011  C  0.0016  v  0.0014  v 

18  6487.0-6601.8  0.0014  C  0.0014  v  0.0019  v 

20  6885.0-6975.7  0.0015  C  0.0016  C  0.0018  C 


TABLE  1.  lntra.sea.son  Dispersion  ami  Interseason  Hanses  in  (b+y)/'2  (conti inns i) 


HD  182572  28  5861.9-5995.6  0.0020  v  0.0019  v  0.0018  C 

21  6218.9-6325.6  0.0014  C  0.0021  V  0.0020  V 

22  6591.9-6718.6  0.0017  C  0.0016  C  0.0018  C 

25  6956.9-7071.6  0.0019  v  0.0016  C  0.0019  v 


[’ABLE  4.  Intraseason  Dispersion  and  Interseason  Ranges  in  (b+y)/ 2  (continued) 


In  addition  to  helping  provide  confirmation  of  intrinsic  variability,  the  ratio 
Sb/sy  is  astrophysically  meaningful:  it  derives  from  the  variation  of  the  disk-averaged 
effective  temperature  of  the  star  caused  by  the  presence  of  cool  dark  spots.  For  stars 
of  solar  temperature,  the  b  filter  lies  shortward  of  the  blackbody  maximum  while 
the  y  filter  is  located  near  the  peak,  so  that  small  temperature  variations  affect  the 
two  filters  unequally:  A b  >  Ay. 

As  a  simple  illustration,  at  6000  K,  a  20  K  decrease  in  disk-averaged  blackbody 
temperature  associated  with  a  spotted  hemisphere  causes  a  0.0185  mag  drop  in  b 
and  an  0.0164  mag  drop  in  y,  yielding  the  ratio  Ab/Ay  =  1.13.  Similarly,  a  20 
K  increase  will  increase  b  and  y  in  the  same  ratio.  Obviously,  spotted  regions  of 
various  temperatures,  spatial  distribution,  and  filling  factors  can  produce  similar 
indistinguishable  effects. 

4.3  Correlations  of  light  curves 

The  principal  statistic  used  to  confirm  variability  is  the  attained  significance 
level  of  the  correlation  coefficient  for  light  curves  containing  one  star  in  common 
(e.g.,  [star  1-star  2]  and  [star  1-star  3]).  In  Table  4, the  designation  V  (variable) 
and  v  (possibly  variable)  indicate,  respectively,  99%  and  95%  significance  of  the 
correlation  coefficient.  Variable  stars  identified  by  this  procedure  are  listed  in  the 
final  column  of  Table  4.  When  it  is  impossible  to  determine  uniquely  from  the  array 
of  correlation  coefficients  (three  for  trios,  six  for  quartets)  which  stars  are  variable, 
as  happens  sometimes  when  more  than  one  star  is  evidently  varying,  the  variability 
designation  is  omitted  even  though  the  Sby  values  of  individual  differential  light 
curves  are  large  enough  that  variability  is  suspected. 

4.4  Frequency  histograms  of  intraseason  variability 

Figure  12  presents  frequency  histograms  of  the  seasonal  standard  deviations, 
.“by,  for  constant  pairs,  possibly  variable  pairs,  and  variable  pairs  (Table  4).  The  data 
divide  evenly  between  variable  and  constant  pairs,  with  only  ~15%  falling  in  the 
possibly  variable  category.  Pairs  with  both  stars  possibly  or  probably  variable  are 
omitted.  Overlap  between  the  histograms  is  small:  the  third  quartile  for  constant 
pairs  occurs  at  0.0026  mag  and  the  first  quartile  for  variable  pairs  is  located  at 
0.0041  mag.  As  a  general  rule,  any  pair  with  Sby  >  0.0025  mag  is  likely  to  include 
a  variable  star,  since  few  of  the  constan1  pairs  have  dispersions  this  large. 

According  to  Figure  12,  half  the  star  pairs  include  a  variable  star;  but  because  of 
redundancy  in  pair  membership,  the  fraction  of  variable  stars  is  actually  fewer  than 
half.  (One  variable  star  in  a  trio  produces  two  entries  in  the  “variable”  histogram 
and  one  in  the  “constant"  histogram.  One  variable  star  in  a  quartet  produces  three 
“variable”  pairs  and  three  "constant”  pairs;  two  variable  stars  in  a  quartet  produce 
five  “variable”  pairs  and  one  “constant”  pair,  etc.)  According  to  the  final  column  of 
Table  4,  35  stars  were  variable  in  two  or  more  seasons,  leading  to  the  conclusion  that 
field  F,  Cl,  and  K  stars  have  a  significant  chance  of  being  found  perceptibly  variable 
on  short  time  scales.  About  half  the  program  stars  and  one-fourth  the  comparison 
stars  were  variable. 
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TABLE  5.  Intraseason  Variability  and  Spectral  Type 


LUMINOSITY 

CLASS 

1  SPECTRAL  TYPE  j 

TOTAL 

Var.  Constant 

Var. 

F 

Constant 

Var. 

G 

Constant 

Var. 

K 

Constant, 

I-III 

1 

2 

3 

9 

5 

2 

9  (41%) 

13 

IVY 

7 

17 

6 

9 

7 

4 

20  (40%) 

30 

Unknown 

1 

11 

1 

6 

3 

3 

5  (25%) 

20 

TOTAL 

9 

30 

10 

24 

15 

9 

34 

63 

(24%) 

(29%) 

(63%) 

(34%) 

Figure  12.  Frequency  histograms  for  star  pairs  exhibiting  definite  variability  at  the 
99%  confidence  level  (top),  at  the  95%  confidence  level  (middle),  or  nonvariability 
( bottom). 
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Intraseason  dispersion  (mag) 


4.5  Consistently  variable  stars 

More  than  one-third  of  the  program  and  comparison  stars  were  found  to  be 
consistently  variable  in  two  or  more  seasons.  Table  5  summarizes  the  number  and 
percentage  distribution  of  variable  stars  by  spectral  type  and  luminosity  classifica¬ 
tion.  About  one-fourth  the  F  and  G  stars  and  two-thirds  the  K  stars  were  variable. 
Surprisingly,  five  of  thirteen  uvby  photometric  standard  stars  are  variable.  Variabil¬ 
ity  at  the  <1%  level  is  therefore  an  extremely  common  phenomenon  among  cool 
stars  regardless  of  luminosity  type. 

Figure  13  shows  the  relation  between  st,y  (averaged  over  all  seasons)  and  B  —  V 
color  for  consistently  variable  program  and  comparison  stars.  The  scale  of  spectral 
types  is  from  Fitzgerald  (1970).  Only  one  star  (HD129333)  fluctuated  by  more  than 
1%.  The  F  and  G  variables  were  predominantly  subgiants  and  dwarfs,  while  the  K 
variables  were  mostly  giants.  Figure  13  includes  seven  uvby  standards,  but  three 
( HD101606,  HD'201091,  and  IID‘201092)  are  marginal  (~0. 002  mag)  detections.  Two 
others  (HD82885,  IID82635),  however,  have  recently  derived  photometric  rotation 
periods  based  on  light  curves  with  amplitudes  greater  than  0.01  mag  (Skiff  and 
Lockwood  1980). 


Figure  13.  Average  intraseason  rms  dispersion,  $byi  for  program  and  comparison 
stars  that  were  variable  in  two  or  more  seasons  as  a  function  of  B  —  V  color.  ( filled 
circles)  -  luminosity  IV  and  V  stars.  ( open  circles)  -  luminosity  class  I,  II,  or  III, 
(stars)  -  luminosity  class  unknown. 
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4.6  Photometrically  derived  rotation  periods 


A  number  of  the  program  stars  have  rotation  periods  determined  earlier  from 
their  UK  variations  (Vaughan  et  al.  1981;  Baliunas  and  Vaughan  1985).  Because 
the  relative  variation  in  the  II K  index,  5,  is  an  order  of  magnitude  greater  than  the 
photometric  variation,  it  is  much  easier  to  measure — typically  A5/5  is  in  the  range 
3-20%,  while  the  rms  b  or  y  variation  in  the  optical  continuum  rarely  exceeds  1%. 

Using  an  algorithm  by  Scargle  (1982),  we  were  nonetheless  able  to  derive  pho¬ 
tometric  rotation  periods  for  nine  stars,  essentially  confirming  the  HK  periods  for 
several,  and  obtaining  original  periods  for  three  others.  Table  6  summarizes  the 
known  HK  and  newly  measured  photometric  rotation  periods.  Theoretical  rotation 
periods,  Pca|c,  derived  from  the  convective  turnover  time  and  B  —  V  color  by  Noyes 
et  al.  (1984)  are  included  for  comparison.  To  maximize  the  signal  to  noise  ratio 
in  this  calculation,  we  first  averaged  together  all  the  differential  magnitudes  having 
the  program  star  in  common;  i.e.,  [(star  1-star  2)  -f  (star  1-star  3)  +  (star  1-star 
4)]/3.  In  Table  6,  the  quantity  Zmax  is  a  measure  of  the  signal- to-noise  ratio  of  the 
photometric  rotation  signal:  if  Zmax  >  10,  it  was  possible  to  derive  photometric 
periods  with  some  confidence. 

5.  INTRINSIC  STELLAR  VARIABILITY  ON  INTERSEASON  TIME  SCALES 

The  detection  of  interseason  variability  is,  in  principle,  straightforward  and 
analogous  to  the  procedure  used  to  analyze  intraseason  variability:  an  analysis  of 
variance  is  used  to  test  the  significance  of  the  differences  between  seasonal  mean 
magnitudes;  and  the  correlation  of  seasonal  mean  light  curves  uniquely  identifies  the 
particular  star  that  is  varying.  However,  this  task  must  be  approached  with  caution 
because,  on  time  scales  of  years,  the  risk  of  systematic  errors  is  ever  increasing. 

Correlations,  or  the  lack  thereof,  among  pairs  of  seasonal  mean  magnitudes 
having  a  subject  star  in  common  were  identified  by  simple  inspection  of  the  light 
curves— the  number  of  annual  data  points  is  presently  too  small  to  justify  formal 
statistics.  Often,  however,  this  method  is  completely  sufficient.  For  example,  in  the 
light  curves  for  the  HD  1835  group  (Figure  4.1),  it  is  obvious  that  the  chromospher- 
ically  active  program  star  (star  1)  is  strongly  variable  from  year  to  year:  the  total 
variation  greatly  exceeds  the  95%  seasonal  confidence  limits  indicated  by  the  error 
bars. 


In  other  cases,  more  formal  statistics  must  prevail.  In  the  HD182572  group 
(Figure  4.25),  for  example,  it  appears  from  inspection  that  star  3  may  be  variable, 
though  its  light  curves  (star  1-star  3,  star  2-star  3)  show  hardly  more  activity  than 
the  constant  pair  (star  1-star  2).  An  analysis  of  variance  confirms  the  suspicion: 
the  light  curves  for  star  1-star  3  and  star  2-star  3  both  indicate  variation  at  a 
confidence  level  >99.9%.  Variation  of  star  1-star  2  is  significant  at  the  99%  level. 

In  cases  of  low-range  variation  hardly  distinguishable  from  noise,  we  must  rely 
upon  significance  tests  in  a  formal  analysis  of  variance  (AOV),  using  the  consis¬ 
tency  among  pairs  of  light  curves  with  a  star  in  common  to  decide  which  star  is 
actually  varying.  Some  groups  present  insoluble  complexity  when  more  than  one 
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TABLE  6.  Summary  of  Rotation  Determinations 


HD 

log  R'tt k 

JD 

2,440,000+ 

n 

sby 

P photorn 

(days) 

^hk  Pc\ic 
(days)  (days)  Zmax 

1835 

-4.42 

6287.0-6432.6 

53 

0.010 

7.7 

7.7 

6.7 

19.5 

25998 

-4.50 

2.6 

35296 

-4.38 

6331.0-6512.6 

63 

0.007 

2.3 

3.3 

10.4 

39587 

-4.40 

6331.0-6512.6 

63 

0.007 

4.9 

5.2 

5.9 

11.4 

82885 

-4.58 

6064.9-6196.7 

54 

0.013 

18.2 

18.1 

18.4 

23.1 

8373.0-6570.7 

55 

0.005 

18.1 

17.5 

20.4 

114710 

-4.73 

6486.9-6570.7 

27 

0.006 

12.4 

12.2 

6.2 

115383 

-4.43 

<0.005 

3.4 

4.9 

115404 

-4.46 

6104.0-6234.7 

39 

0.006 

11.0: 

18.8 

15.5 

7.2 

129333 

-4.23 

6104.0-6250.7 

53 

0.021 

3.0 

17.2 

6487.0-6639.7 

52 

0.014 

2.7 

11.7 

131156  AB 

-4.36 

6104.0-6250.7 

43 

0.006 

6.1 

6.2b 

7.7 

12.5 

6487.0-6601.7 

33 

0.007 

6.4 

12.6 

6882.0-6979.7 

35 

0.008 

6.2 

9.2 

149661 

-4.54 

0.006 

21.3 

17.7 

160346 

-4.71 

<0.005 

33.5 

37.0 

176095 

-4.68 

31 

0.009 

2.8 

10.4 

190007 

-4.59 

5861.9-6014.6 

6 

0.008 

6218.9-6375.6 

37 

0.004 

35.0 

29.3 

25.2 

11.6 

6587.9-6731.6 

48 

0.006 

28.0 

9.7 

5861.9-6998.4 

147c 

0.008 

26.3 

14.8 

201091 

-4.80 

0.003 

37.9 

40.6 

9.1 

201092 

-4.91 

48.0 

50.0 

NOTES: 

*  Noyes  et  al.  1981. 
b6.2-6.7  (Baliunas  et  al.  1985) 
cFour  years  combined. 
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star  is  evidently  varying;  e.g.,  IID161239  (Figure  4.23),  HD201091/HD201092  (Fig¬ 
ure  4.28). 

On  the  bottom  row  of  data  for  each  group.  Table  4  summarizes  the  results 
of  the  formal  statistical  tests  for  interseason  variability.  The  peak-to-peak  range 
(b  and  y  averaged)  for  each  star  pair  is  listed  along  with  the  designation  V,  v,  or 
C  according  to  the  outcome  of  the  AOV  for  that  pair:  V  means  that  the  attained 
significance  (via  the  F  statistic)  for  variation  in  the  annual  mean  values  is  99%  or 
greater;  v  indicates  a  significance  greater  than  95%,  and  C  indicates  no  significant 
variability.  By  inspecting  the  decision  for  each  pair,  the  identity  of  the  variable 
star(s)  is  usually  evident  and  is  listed  in  the  final  columnm  of  Table  4.  A 
denotes  variation  in  excess  of  0.005  mag  (0.5%)  peak-to-peak —these  are  the  most 
certain  detections. 

Figure  14  reports  the  outcome  of  the  analysis  of  variance:  the  upper  histogram 
pertains  to  pairs  of  stars  that  show  evidence  of  variability  at  the  99%  confidence  level, 
the  center  histogram  includes  pairs  in  the  95-99%  range,  and  the  bottom  histogram 
includes  nonvariable  pairs.  According  to  this  procedure,  over  half  the  pairs  show 
evidence  of  some  variability.  As  there  remains  some  possibility  of  systematic  error 
contributing  to  the  observed  variability,  this  statistic  alone  is  insufficient;  and  we 
require  further  evidence  before  deciding  that  a  particular  star  is  variable,  including 
correlation  of  the  appropriate  pairs  of  light  curves,  and  the  existence  of  a  pair  of 
stars  in  each  group  deemed  “constant,”  regardless  of  the  AOV. 

In  all  groups,  but  especially  in  quartets,  one  pair  of  stars  has  to  be  chosen 
as  “constant,”  regardless  of  the  outcome  of  the  formal  AOV,  to  provide  a  basis 
for  assessing  the  variations  of  the  other  stars.  Sometimes  this  pair  itself  is  suspect 
because  the  total  amplitude  of  variation  is  clearly  larger  than  that  seen  in  some  of 
the  best-behaved  groups.  Most  likely,  this  simply  indicates  that  no  pair  of  stars  in 
that  group  is  truly  constant— or  it  could  mean  that  a  systematic  error  of  unknown 
origin  is  contaminating  the  data  for  that  particular  group.  We  consider  this  problem 
further  below. 

Analogous  to  Figure  13  for  intraseason  variations,  we  show  in  Figure  15  the 
relation  between  peak-to-peak  three-  or  four-year  range  of  seasonal  mean  magnitudes 
( b  and  y  averaged)  and  B  -  V  color.  The  median  range  is  0.008  mag,  and  there 
are  no  trends  with  luminosity  or  spectral  type.  Ranges  at  the  0.002  mag  level  are 
marginally  significant  but  may  ultimately  be  confirmed  by  further  observations. 

5.1.  The  least-variable  stars 

Having  established  (Table  5)  that  one-third  of  the  entire  sample  and  two-thirds 
of  the  K  stars  are  consistently  variable  on  intraseason  time  scales,  we  now  consider 
the  statistics  of  long-term  variability.  The  least-variable  stars  are  provisionally  so 
designated  and  listed  by  differential  pairs  in  Table  7  if  they  (1)  show  a  total  range  of 
the  averaged  b  and  y  seasonal  mean  magnitudes  <0.003  mag  and  (2)  have  an  average 
intraseason  rms  dispersion  <0.003  mag.  Twenty-six  pairs  of  stars,  including  40 
individual  stars  (12  appear  in  two  or  more  pairs)  qualify  for  this  list.  Even  with  this 
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Figure  14.  Frequency  histograms  of  peak-to-peak  range  of  seasonal  mean  mag¬ 
nitudes  for  pairs  of  stars  judged  variable  at  the  99%  confidence  level  ( top  panel), 
possibly  variable  at  the  95-99%  confidence  level  ( middle  panel),  or  nonvariable  ( bot¬ 
tom  panel). 

severe  restriction  upon  the  total  range  of  variability,  17  pairs  of  stars  appear  to  be 
significantly  (>95%  probable)  variable  according  to  the  formal  analysis  of  variance, 
leaving  only  9  pairs  showing  no  evidence  whatsoever  of  long-term  variability.  These 
18  “constant”  stars  are  predominantly  F  stars;  and  the  sole  MK-classified  giant  in 
the  list,  HD73596  (F5III),  may  fluctuate  slightly  each  year.  Another  likely  giant, 
HD140716  (gG9),  the  coolest  of  the  “constant”  stars,  is  also  suspected  of  varying 
slightly  each  year. 
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Interseason  range  (mag) 


Figure  15.  Iuterseason  range  for  program  and  comparison  stars  that  were  judged 
variable  according  to  the  analysis  of  variance  and  correlated  light  curves,  (filled 
circles )  -  luminosity  class  IV  and  V  stars,  (open  circles)  -  luminosity  class  I,  II,  or 
III,  (stars)  -  luminosity  class  unknown. 

5.2  Possible  long-term  systematic  effects 

The  most  likely  long-term  systematic  error  source  we  could  imagine  has  to  do 
with  the  slowly  changing  color  response  of  the  photometer  that  might  accompany 
the  aging  of  the  photomultiplier  or  the  interference  filters.  Such  an  effect,  if  present, 
would  be  most  easily  seen  as  a  drift  in  the  differential  magnitudes  of  two  stars  widely 
different  in  B  —  V  color  (cf.  Figure  8).  The  group  of  least-variable  stars  (Table  7) 
contains  pairs  of  stars  whose  B  -  V  colors  differ  by  up  to  0.6  mag.  On  Figure  16 
are  plotted  the  colors  of  the  hotter  members  of  the  various  pairs  as  a  function  of 
color  difference,  using  symbols  whose  size  represents  the  range  of  variation  (larger 
symbols  for  larger  range,  up  to  0.003  mag).  There  is  no  tendency  for  the  star  pairs 
with  the  larger  ranges  of  variability  to  cluster  anywhere  on  this  diagram;  therefore 
we  conclude  that  systematic  color-related  effects  are  absent. 

5.3.  The  large-range  variable  stars:  incipient  cycles? 

Stars  exhibiting  “large”  variations  (i.e.,  >0.5%  peak-to-peak  in  the  light  curve 
of  annual  mean  magnitudes)  over  the  four-year  course  of  this  program  (cf.  Table  4) 
comprise  about  one-fourth  of  the  total  sample  and  are  found  predominantly  among 
the  program  stars.  These  stars  are  listed  in  Table  8  along  with  an  indication  of  their 
intraseason  variability:  a  “V”  signifies  consistent  variability  in  two  or  more  seasons, 
and  a  “C”  signifies  nonvariability  in  most  seasons.  Typical  values  of  the  intraseason 
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TABLE  7.  The  Least- Variable  Stars 


HD 

Spectral 

Type 

Range 

(mag) 

AOV 

Average 

s.d. 

Note 

2488 

1388 

F5 

GO  V 

0.0013 

C 

0.0027 

10697 

1 1326 

G5  IV 

G5 

0.0009 

C 

0.0022 

13121 

13683 

GO  V 

F0 

0.0020 

C 

0.0021 

uvby,  UK  std,  log  =  -5.21 

18256 

18404 

F6  V 

F5  IV 

0.0011 

C 

0.0017 

lo§  R'hk  =  ~4-79 

18404 

17659 

F5  IV 

F8 

0.0024 

V 

0.0019 

Variable? 

23256 

22679 

F2 

G5 

0.0019 

V 

0.0020 

Three  seasons  only 

76572 

73596 

F6  V 

F5  III 

0.0024 

V 

0.0020 

HK  std,  log  =  —4.94,  variable 

Var?  each  year 

76572 

78234 

F6  V 

F2  V 

0.0020 

V 

0.0019 

7.3596 

78234 

F5  III 

F2  V 

0.0007 

c 

0.0022 

Var?  each  year 

83951 

83525 

F3  V 

F5 

0.0018 

V 

0.0019 

Three  seasons  only 

103095 

103520 

G8  VI 

K0  III 

0.0015 

V 

0.0018 

uvby  std,  log  R'HK  =  —4.87 

Var?  in  two  seasons 

103095 

101606 

G8  VI 

F4  V 

0.0028 

V 

0.0016 

uvby  std,  variable 

124570 

125451 

F6  IV 

F5  IV 

0.0025 

V 

0.0017 

HK  std,  log  R'hk  =  —5.14 

129390 

131330 

F2 

F8 

0.0014 

c 

0.0027 

dispersion  (b  and  y  averaged),  and  the  interseason  observed  range  are  listed. 

The  sixteen  program  stars  listed  in  Table  8  are  for  the  most  part  consistently 
variable  on  intraseason  time  scales  as  well.  Only  four  stars,  those  with  the  smallest 
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TABLE  7.  The  Least- Variable  Stars  continued. 


HD 

Spectral 

Type 

Range 

(mag) 

AOV 

Average 

s.d. 

Note 

1  13701 
112091 

GO  +Va 
K1  Va 

0.0022 

V 

0.0016 

uvby,  MK,  II K  std,  log  =  -5.02 

Variable 

113761 

140716 

GO  -f  Va 
gC9 

0.0014 

c 

0.0020 

Slightly  variable? 

142091 

140716 

Kl  Va 
gC9 

0.0029 

V 

0.0017 

Variable 

157856 

157347 

F3  V 

G5  IV 

0.0025 

V 

0.0025 

log  =  -4.70,  variable? 

156635 

157347 

F8 

G5  IV 

0.0015 

c 

0.0022 

176095 

175515 

F5  IV 

K0  III 

0.0025 

V 

0.0023 

log  R'hk  =  -4.68 

182572 

180868 

G8  IV 

F0  IV 

0.0017 

V 

0.0018 

MK  std,  log  =  -5.03 

200031 

201092 

G5  III 

K7  V 

0.0015 

V 

0.0020 

uvby,  MK  std,  log  =  -4.91 

200031 

200577 

G5  III 

G8  III 

0.0027 

V 

0.0019 

201092 

200577 

K7  V 

G8  III 

0.0021 

V 

0.0020 

215704 

216175 

K0 

G5 

0.0026 

V 

0.0023 

<  5  >=  0.25. 

216048 

217232 

F0 IV- V 
A8  V 

0.0014 

c 

0.0023 

Both  slightly  variable? 

(and  most  Sun-like)  log  R'HK  values,  failed  to  exhibit  consistent  intraseason  vari¬ 
ability.  The  eight  comparison  stars  were  evenly  divided,  with  just  half  consistently 
varying  within  two  or  more  seasons.  Because  the  sample  here  is  small,  it  seems 
imprudent  to  attempt  to  make  a  distinction  between  the  statistics  of  the  program 
and  comparison  stars  in  Table  8,  even  though  the  latter,  by  and  large,  is  comprised 
of  more-luminous  stars. 

Do  the  light  curves  of  these  stars  contain  clues  to  the  nature  of  the  variability? 
Can  random  variability  be  distinguished  from  segments  of  cycles?  To  consider  these 
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(B-V)  color  of  the  hotter  star 


Figure  16.  B-V  colors  of  the  hotter  members  of  the  pairs  of  least- variable  stars 
from  Table  7  as  a  function  of  color  difference.  The  size  of  the  symbols  is  related  to 
the  observed  range  of  interseason  variability,  in  the  range  0.0007  mag-0.0029  mag. 

questions,  we  examined  the  statistics  of  artificial  “light  curves ”  consisting  of  sets 
of  four-point  time  series  generated  from  normally  distributed,  consecutive  random 
numbers.  Fourteen  percent  of  the  artificial  “light  curves”  thus  generated  were  mono¬ 
tonic  (steadily  rising  or  falling),  52%  contained  a  single  maximum  or  minimum,  and 
34%  contained  two  extrema — one  maximum  and  one  minimum.  Thus  66%  of  the 
random-noise  light  curves  had  characteristics  consistent  with  cycle  segments  where 
the  periods  are  substantially  longer  than  four  years. 

However,  turning  to  the  real  data,  only  13%  of  the  observed  light  curves  (Table 
8)  contained  two  extremes,  23%  were  monotonic,  and  64%  contained  just  a  single 
maximum  or  minimum.  Thus  86%  were  consistent  with  undersampled  cyclic  varia¬ 
tion,  compared  with  06%  of  the  artificial  light  curves,  and  the  fraction  of  observed 
double-extrema  light  curves  was  less  than  half  the  fraction  generated  artificially. 
These  statistics  suggest  that  the  observed  light  curves  are  slightly  more  consistent 
with  short  segments  of  some  kind  of  cyclic  phenomenon  than  with  random  noise  on 
a  one-year  time  scale.  This  conclusion  is,  of  course,  highly  preliminary  because  of 
the  short  time  base  of  available  data;  and  it  is  prejudiced  by  the  expectation  that 
luminosity  cycles,  if  detected,  will  be  commensurate  in  length  with  the  magnetic 
activity  cycles. 
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TABLE  8.  Strongly  Varying  Stars 


HD 

Spectral 

Type 

Intra 

Log  ft'HK  Season 

Typ.  rms 
Dispersion 
(mag) 

Peak  Peak 
Range 
(mag) 

Number 

of 

Extremes 

Program  Stars 

1835 

G2V 

-4.42 

V 

0.009 

0.013 

2 

35296 

F8V 

-4.38 

V 

0.007 

0.010 

1 

39587 

GOV 

-4.44 

V 

0.006 

0.015 

1 

75332 

F7Vn 

-4.44 

V 

0.007 

0.018 

82885 

G8IV-V 

-4.58 

V 

0.007 

0.008 

1 

82635 

G8.5III 

V 

0.005 

0.005 

1 

115383 

GOVs 

-4,13 

V(?) 

0.003 

0.015 

1 

115404 

K1V+M1V 

-4.46 

V 

0.007 

0.013 

1 

120136 

F6IV 

-4.76 

c 

0.003 

0.011 

0 

129333 

GO 

-4.23 

V 

0.017 

0.048 

0 

149661 

K2V 

-4.54 

V 

0.006 

0.010 

1 

152391 

G7V 

-4.45 

V 

0.007 

0.022 

1 

158614 

G9IV-V 

-4.99 

c 

0.002 

0.006 

2 

161239 

G2IIIb 

-5.14 

c 

0.004 

0.005 

1 

185144 

K0V 

-4.75 

c 

0.002 

0.007 

? 

190007 

K4 

-4.59 

V 

0.007 

0.010 

1 

Comparison  Stars 

13611 

G6II-III 

c 

0.002 

0.017 

0 

112989 

G91II 

V 

0.006 

0.013 

1 

126126 

F8V+G1V 

V 

0.004 

0.018 

1 

129972 

G8.5III 

c 

0.007 

0.011 

0 

131511 

K2V 

V 

0.003 

? 

0 

160823 

GO 

c 

0.004 

0.005 

2 

160935 

F5 

c 

0.003 

0.007: 

1 

215427 

K5 

V 

0.010 

0.032 

1 

G.  PHOTOMETRIC  VARIABILITY  AND  CHROMOSPHERIC’  ACTIVITY 

Figure  17  shows  the  seasonal  values  of  the  dispersion,  st,y,  for  each  of  the  pro¬ 
gram  stars  as  a  function  of  the  chromospheric  emission  ratio  log  One  symbol 

is  plotted  for  each  season’s  observations.  Generally,  the  magnitude  of  var‘ ability 
each  season  is  comparable,  leading  to  the  observed  clustering  of  points  on  the  dia¬ 
gram,  and  the  conclusion  that  the  amount  of  spot  coverage  may  change  only  slowly 
with  time.  The  G8IV-V  program  star  HD82885  provides  an  example  where  the 
?by  values  have  showed  a  clear  evolution:  in  four  consecutive  seasons,  S{,y  =  0.0066, 
0.0127  (with  a  rotational  period  determined),  0.0057,  and  0.0035  mag.  Obviously, 
it  is  important  to  determine  if  the  changes  in  seasonal  activity  are  in  any  way  cyclic 
and  if  the  intraseason  variability  can  be  related  to  the  average  brightness. 
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Chromospheric  emission  ratio  log  R'^ 


Figure  17.  Intraseason  standard  deviation,  Sby,  as  a  function  of  the  chromospheric 
emission  ratio  log  R'HK.  Each  star  is  plotted  once  for  each  season  in  which  a  vari¬ 
ability  decision  could  be  made.  ( open  circles)  -  constant,  (star)  -  possibly  variable, 
(filled  circles)  -  variable. 

In  Figure  17,  the  onset  of  consistent  intraseason  activity  occurs  rather  suddenly 
at  a  value  of  log  ft[(K  ~  -4.6,  about  twice  (in  linear  units)  the  solar  value  of  -4.94. 
In  this  regard,  the  data  seem  rather  significantly  grouped:  there  is  only  one  program 
star  with  log  R'HK  >  -4.6  with  an  observed  intraseason  dispersion  less  than  0.003 
mag.  Below  -4.6,  the  regime  of  the  true  solar  analogs,  the  opposite  holds:  no  star 
is  more  variable  than  about  0.004  mag.  We  note  that  the  Sun  itself,  if  observed  in 
this  manner,  would  be  seen  as  quiescent,  with  a  dispersion  well  under  0.001  mag. 

Figure  18  shows  the  three-  or  four-year  range  of  variation  for  the  program 
stars  as  a  function  of  the  chromospheric  emission  ratio,  log  R'HK.  The  error  bar 
tails  indicate  the  95%  confidence  interval:  if  the  error  bar  lies  fully  above  zero,  the 
star  is  considered  variable.  By  this  criterion,  most  of  the  program  stars  seem  to  be 
variable.  Note  the  position  of  the  Sun  on  the  diagram,  based  on  an  estimated  cyclic 
amplitude  of  0.1%. 

Unlike  intraseason  variabiliy,  the  ranges  of  interseason  variation  do  not  show  a 
distinct  grouping,  except  for  the  general  rise  with  increasing  log  R'HK-  In  particular, 
among  the  “active”  stars  (log  /?jiK  >  — 1.6),  some  stars  maintain  a  very  small  range 
of  observed  long-term  variation. 
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Peak-to-peak  range  of  st.  >nal  mean  mags 


Figure  18.  Interseason  range  of  variation  over  four  years  plotted  as  a  function  of 
the  chromospheric  emission  ratio  log  R'HK .  The  error  bar  “tails”  indicate  the  95% 
confidence  interval  for  the  range:  if  the  bars  lie  above  zero,  the  star  is  considered 
variable. 


Intraseason  and  interseason  variability  are  intimately  related,  at  least  among 
the  program  stars.  Of  the  16  program  stars  considered  “strongly  variable”  because 
their  four- year  ranges  exceed  0.005  mag,  all  but  four  have  values  of  log  /?hK  >  —4.6, 
and  it  is  these  four  alone,  those  most  like  the  Sun  in  terms  of  emission  ratio,  that  fail 
to  display  intraseason  variability  as  well.  These  four  stars  are  HD120136  (F6IV), 
HD158614  (G9IV-V),  HD161239  (G2IIIb— the  sole  giant),  and  HD190007  (K4). 

7.  SUMMARY  AND  CONCLUSIONS 

To  characterize  the  luminosity  variations  of  stars  similar  to  the  Sun  but  cov¬ 
er:  ng  a  range  of  stellar  age,  we  have  obtained  over  2,000  nightly  observations  of 
36  program  stars  during  the  interval  March  1984  through  December  1987.  Chro¬ 
mospheric  activity  in  the  program  stars  was  previously  monitored  for  a  decade  by 
0.  C.  Wilson,  who  discovered  that  many  of  the  presumably  old  solar  analog  stars 
exhibited  cycles  of  magnetic  activity  similar  to  that  of  the  Sun.  The  program  stars 
were  selected  to  bracket  the  Sun  in  several  observable  quantities:  (1)  temperature, 
as  suggested  by  B  -  V  color  and  MK  spectral  type;  (2)  age,  as  suggested  by  atomic 
abundances  and  chromospheric  activity;  (3)  rotation  rate,  as  indicated  by  periodic¬ 
ities  in  the  strength  of  chromospheric  emission;  (4)  mass,  as  inferred  from  a  variety 
of  generally  accepted  astrophysical  relationships;  and  (5)  known  cyclic  behavior  of 
magnetic  activity. 
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The  stars  were  measured  differentially  at  472  and  551  mil  with  respect  to 
two  or  three  nearby  comparison  stars  of  similar  brightness  and  color,  with  the  aim 
of  obtaining  the  most  precise  observations  possible.  Accordingly,  we  have  devoted 
considerable  attention  to  "quality  control":  all  observations  were  made  using  a  stan¬ 
dard  time-tested  procedure  that  has  built-in  checks  of  consistency  and  precision  and 
which  produces  data  sets  of  equal  statistical  weight,  unbiased  with  regard  to  a  priori 
assumptions  about  which  star(s)  within  groups  may  be  variable.  Though  not  ob¬ 
tained  using  a  fully  randomized  design,  the  observations  seem  to  be  free  of  temporal 
artifacts  and  the  results  demonstrate  a  precision  that  exceeds  the  expectations  of 
the  original  experimental  objective. 

The  principal  conclusions  of  this  study  are: 

1.  About  one-half  the  program  stars  and  one-IA..  h  the  comparison  stars  are 
demonstrably  variable  on  intraseason  time  scales,  as  it,  ..  .tied  by  («)  the  magnitude 
of  the  night-to-night  dispersion  of  differential  magnitudes,  (6)  the  correlation  of 
pairs  and  trios  of  light  curves  having  a  variable  star  in  common,  (c)  the  repeated 
occurrence  of  large  dispersions  in  two  or  more  independent  observing  seasons,  and 
(d)  the  ratio  of  dispersions  in  the  two  filter  passbands,  consistent  with  a  “starspot” 
model  of  irradiance  variation.  The  boundary  between  variable  and  nonvariable  stars 
is  about  0.0(J3  mag  (0.3%)  rms  in  the  cycle- to- cycle  dispersion  of  averaged  b  and 
y  differential  magnitudes  averaged  together.  The  range  of  recorded  variability  is 
from  0.002-0.02  mag,  with  a  median  of  0.000  mag  (compared  with  a  median  value 
of  0.002  mag  for  constant  stars). 

2.  The  incidence  of  intraseason  variability  is  essentially  independent  of  lumi¬ 
nosity  class — about  40%  for  classes  I-III  and  IV-V  alike — but  is  about  twice  as  likely 
among  K  stars  (63%)  as  among  F  and  G  stars,  regardless  of  luminosity.  However 
the  amplitude  of  observed  variation  among  the  consistently  mumMo  stars  is  essen¬ 
tially  independent  of  B  -  V  color,  even  though  it  may  arise  from  different  physical 
phenomena,  depending  upon  temperature  and  luminosity. 

3.  The  determination  of  rotation  periods  is  relatively  more  difficult  from  pho¬ 
tometry  than  from  the  chromospheric  UK  index  because  the  signal  amplitudes  are 
much  smaller.  Nonetheless,  rotation  periods  were  determined  for  nine  stars,  includ¬ 
ing  three  for  which  no  UK  rotation  rates  were  available.  For  four  stars,  rotation 
periods  were  measured  in  at  least  two  seasons,  but  the  data  are  insufficient  to  tell 
if  the  season- to-season  differences  are  significant. 

4.  Using  principally  an  analysis  of  variance,  most  of  the  differential  light  curves 
exhibited  statistically  significant  interseason  variations  with  three-  or  four-year  peak- 
to-peak  ranges  from  0.002  0.05  mag  (0.2-5%).  With  only  four  data  points  per 
light  curve,  a  formal  correlation  analysis  is  premature;  however,  in  many  cases, 
simple  inspection  of  the  ensemble  of  light  curves  is  sufficient  to  identify  variable 
stars  unambiguously.  Twenty-four  stars  (including  16  of  the  36  program  stars) 
exhibited  large  enough  ranges  (>0.005  mag)  in  two  or  more  light  curves  to  indicate 
variability  beyond  reasonable  doubt.  Such  large-range  variability  was  relatively 
uncommon  among  the  comparison  stars,  occurring  in  only  8  out  of  a  total  of  61 
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stars  observed  in  three  or  more  seasons.  By  usual  standards,  therefore,  a  great 
majority  of  comparison  stars  are  stable  enough  to  be  used  as  standard  stars  for 
conventional  photometry,  even  though  they  may  be  demonstrably  variable  by  small 
amounts.  Nearly  all  the  stars  that  vary  strongly  on  interseason  time  scales  are 
consistently  variable  on  intraseason  time  scales,  with  the  notable  exception  of  four 
program  stars  with  near-solar  values  of  log  (see  below). 

5.  Variability  of  the  program  stars  is  a  strong  function  of  the  chromospheric 
emission  ratio  A’JiK'  and  it  tends  to  repeat  at  about  the  same  level  from  year  to 
year.  Consistent  intraseason  variability  is  observed  in  stars  with  emission  ratios 
about  twice  the  solar  value  and  greater  (log  R\lK  >  -1.6).  luterseason  variability 
likewise  is  a  strong  function  of  emission  ratio,  but  four  stars  with  near-solar  values  of 
log  have  shown  total  ranges  greater  than  0.005  mag,  even  though  intraseason 
variability  was  absent. 

6.  The  identification  of  cyclic  behavior  is  not  yet  possible  for  program  stars 
because  the  time  base  is  insufficient.  However,  it  is  significant  that  of  the  16  program 
stars  exhibiting  total  four-year  ranges  greater  than  0.005  mag,  twelve  of  the  light 
curves  are  consistent  with  cyclic  variation:  either  they  rise  or  fall  monotonically,  or 
they  exhibit  at  most  one  extremum — i.e.,  “flickering’’  is  rarely  seen.  This  strongly 
suggests  that  further  observation  will  reveal  cyclic  variation. 
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APPENDIX 


TYPICAL  RAW  DATA  FILE  AS  RECORDED  AT  THE  TELESCOPE 

Entries  to  the  right  of  the  colons  in  the  eight-line  header  are  system  defaults 
and  information  provided  by  the  observer.  Additionally,  remarks  may  be  inserted 
into  the  data  stream  at  any  point  and  are  keyed  by  the  time  of  entry  (this  file 
happens  to  have  none). 

One  data  line  is  produced  for  each  10-second  integration.  The  columns,  re¬ 
spectively,  are  UT  month,  day,  year,  hour,  minute,  second,  star  identification  (in 
this  case  2,  1,  4,  or  3),  filter  number  (2  =  Sr90  source,  3  =  y,  4  =  6),  whether 
the  observation  is  a  “star”  or  a  “sky”  (ST  or  SK),  a  numerical  code  for  the  same 
information,  and,  finally,  the  number  of  counts  recorded. 

The  order  of  integrations  is  the  same  as  indicated  in  §2.3  (page  17):  Dark, 
std  source,  std  source,  sky,  sky,  star,  star,  star  ....  Note  at  04:51:51  the  filter  is 
changed  from  y  to  6.  and  at  05:06:10  from  b  back  to  y.  The  total  elapsed  time 
for  four  cycles,  y,h,6,y  is  31  minutes,  yielding  a  data  collection  efficiency  of  75% 
(a  remarkable  feat  for  a  manually  slewed  telescope!). 
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OBSERVER:  BAS  ON  LVS099,  870328.Q02 
I’ HOG  RAM:  LVSSS.HD  82885 

WEATHER:  CLEAR  >  WEAK  TROUGH,  CALM,  DARK 

COMMENTS:  DIAPHRAGM  4,  SEEING  ~8" 

TEMPERATURE:  24F,  BOX  @  69F 

TELESCOPE:  21  INCH 

INTEG.  TIME:I0  SEC 

S/N:  300  NPAUSS:  10  NPAUSB:  10 
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03  28  1987  05  03  34  3 

03  28  1987  05  03  44  3 

03  28  1987  05  03  56  3 

03  28  1987  05  04  06  3 

03  28  1987  05  04  23  3 

03  28  1987  05  04  33  3 

03  28  1987  05  04  43  3 

03  28  1987  05  05  13  2 

03  28  1987  05  05  23  2 

03  28  1987  05  05  33  2 

03  28  1987  05  05  46  2 

03  28  1987  05  05  56  2 

03  28  1987  05  06  10  2 

03  28  1987  05  06  20  2 

03  28  1987  05  06  39  2 

03  28  1987  05  06  49  2 

03  28  1987  05  06  59  2 

03  28  1987  05  07  24  1 

03  28  1987  05  07  35  i 

03  28  1987  05  07  45  1 

03  28  1987  05  07  57  1 

03  28  1987  05  08  07  1 

03  28  1987  05  08  26  1 

03  28  1987  05  08  36  1 

03  28  1987  05  08  46  1 


4  ST  1  1043625 

4  ST  1  1050818 

4  ST  1  1037219 

4  SK  2  1361 

4  SK  2  1357 

4  ST  1  1045662 

4  ST  1  1041850 

4  ST  1  1042583 

4  ST  1  297718 

4  ST  1  296424 

4  ST  1  297957 

4  SK  2  1305 

4  SK  2  1262 

4  ST  1  296825 

4  ST  1  295665 

4  ST  1  295704 

4  ST  1  647121 

4  ST  1  645897 

4  ST  1  643341 

4  SK  2  1417 

4  SK  2  1388 

4  ST  1  646031 

4  ST  1  646196 

4  ST  1  644847 

4  ST  1  2102190 

4  ST  1  2113119 

4  ST  1  2112761 

4  SK  2  1485 

4  SK  2  1396 

3  SK  2  2065 

3  SK  2  2047 

3  ST  1  1699506 

3  ST  1  1703262 

3  ST  1  1704188 

3  ST  1  786593 

3  ST  1  790224 

3  ST  1  783932 

3  SK  2  2028 

3  SK  2  2000 

3  ST  1  782534 

3  ST  1  782951 

3  ST  1  782128 


03  28  1987 

05  09  13 

4 

03  28  1987 

05  09  23 

4 

03  28  1987 

05  09  33 

4 

03  28  1987 

05  09  45 

4 

03  28  1987 

05  09  55 

4 

03  28  1987 

05  10  11 

4 

03  28  1987 

05  10  21 

4 

03  28  1987 

05  10  31 

4 

03  28  1987 

05  11  02 

3 

03  28  1987 

05  11  12 

3 

03  28  1987 

05  1 1  22 

3 

03  28  1987 

05  11  34 

3 

03  28  1987 

05  11  44 

3 

03  28  1987 

05  12  01 

3 

03  28  1987 

05  12  11 

3 

03  28  1987 

05  12  22 

3 

03  28  1987 

05  12  51 

2 

03  28  1987 

05  13  01 

2 

03  28  1987 

05  13  11 

2 

03  28  1987 

05  13  23 

2 

03  28  1987 

05  13  33 

2 

03  28  1987 

05  13  47 

2 

03  28  1987 

05  13  57 

2 

03  28  1987 

05  14  13 

2 

3 

ST 

1 

197812 

3 

ST 

1 

198589 

3 

ST 

1 

197407 

3 

SK 

2 

1971 

3 

SK 

2 

1923 

3 

ST 

1 

197126 

3 

ST 

1 

198257 

3 

ST 

1 

195730 

3 

ST 

1 

400082 

3 

ST 

1 

399779 

3 

ST 

1 

403006 

3 

SK 

2 

2093 

3 

SK 

2 

2173 

3 

ST 

1 

401794 

3 

ST 

1 

402532 

3 

ST 

1 

400210 

3 

ST 

1 

1700011 

3 

ST 

1 

1701801 

3 

ST 

1 

1700465 

3 

SK 

2 

2092 

3 

SK 

2 

2158 

2 

ST 

1 

389807 

2 

ST 

1 

391184 

2 

SK 

2 

166 

REDUCED  DATA  OUTPUT,  SAMPLE  PRINTED  FILE 

On  the  first  page,  header  information  from  the  raw  data  file  is  repeated, 
amended  to  include  input  and  output  file  names,  the  name  of  the  file  in  which 
star  positions  are  found,  assumed  value  of  extinction,  and  assumed  positions  of 
the  stars  in  the  group.  The  tabulation  following  gives  the  fractional  time  of  day. 
an  internal  star  identification  number,  filter  number  (S  =  standard  source,  4  = 
y,  2  =  6),  instrumental  magnitude  !•  u  each  sequence  of  three  10-second  inte¬ 
grations  (corrected  for  sky  background,  pulse  pair  dead  time,  and  atmospheric 
extinction),  and  in  the  final  column  an  indication  of  outlying  data  points  (if  any). 

The  second  page  includes  in  the  upper  section  a  tabulation  of  the  differential 
magnitudes  and  differential  air  masses  for  each  pairwise  combination  and  for 
each  cycle  of  observation.  These  values  are  reported,  cycle  by  cycle,  to  the  stored 
output  file  named  in  the  header  on  the  previous  page,  and  are  used  for  further 
analysis. 

Finally,  a  summary  at  the  bottom  of  the  page  provides  averaged  differential 
magnitudes  for  that  particular  data  set. 


t 
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PHOTOMETRY  OF  HD  82885 

RAW  DATA  FILE:  870328.Q02  ’’OSITION  FILE:  LVSSS.POS 
INPUT  FILE:  870328.S04  (  rOUP  IDENT.=  4 

OUTPUT  FILE:  870328.T04  EXTINCTION:  KB  =  0.230  I\'V=  0.360  KY=  0.170 

OBSERVER:  BAS  ON  I.VS099,  870328. Q02 
PROGRAM:  LYSSS.HD  82885 

WEATHER:  (  ,  ZAR  >  WEAK  TROUGH,  CALM,  DARK 
COMMENTS:  DIAPHRAGM  4,  SEEING  -8" 

I  E.M  I'ERATURE:  24 F,  BOX  @  ODE 
TELESCOPE:  21  INCH 
INTLG.  TIME:  10  SEC 

S/N:300  NPAUSS:10  NPAUSB:10 
COORDINATES  OF  OBJECT  1  9  35.7  35  49 

COORDINATES  01  OBJECT  2  9  3  E2  36  24 

COORDINATES  Ol  OBJECT  3  9  42.7  35  6 

COORDINATES  OI  OBJECT  4  9  40.0  35  20 


1987  3  28 


AC.  DAY 

♦ 

F 

D 

A.M. 

S 

- 

6.0204 

0.1974 

2 

4 

4.2140 

1.0003 

(i  1980 

1 

4 

5.0930 

1.0001 

0  1987 

1 

4 

5.0936 

1.0002 

0.1992 

4 

4 

6.6005 

1.0000 

0.2000 

4 

4 

6.6000 

1.0000 

0.2005 

3 

4 

5,8288 

1.0000 

0.2012 

3 

4 

5.8300 

1.0000 

0.2017 

2 

4 

1.2155 

1.0008 

0.2028 

2 

2 

3.9478 

1.0010 

0.2033 

1 

2 

4.7218 

1.0008 

0.2040 

1 

2 

4.7225 

1.0009 

0.2046 

4 

2 

6.0902 

1.0004 

0.2053 

4 

2 

6.0930 

1.0005 

0.2058 

3 

2 

5.2414 

1.0003 

0.2065 

3 

2 

5.2436 

1  0004 

0.2071 

2 

2 

3.9463 

1 .0022 

STAR- SKY 

N  COUNT  >  2  SIGMA  FROM  MEAN 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


% 


-  99  - 


0.2078 

2 

2 

3.9502 

1.0024 

3 

0.208-1 

1 

2 

4.7186 

1.0021 

3 

0.2090 

1 

2 

4.7191 

1.0023 

3 

0.2096 

4 

2 

6.0895 

1.0015 

3 

0.2103 

4 

2 

6.0942 

1.0017 

3 

0.2108 

3 

2 

5.2439 

1.0014 

3 

0.2115 

3 

2 

5.2434 

1.0015 

3 

0.2121 

2 

2 

3.9479 

1.0041 

3 

0.2131 

2 

4 

4.2436 

1.0046 

3 

0.2136 

1 

4 

5.0880 

1.0042 

3 

0.2113 

1 

4 

5.0940 

1.0045 

o 

0.2148 

4 

4 

6.5978 

1.0034 

3 

0.2155 

4 

4 

6.6027 

1.0037 

3 

0.2161 

3 

4 

5.8253 

1.0031 

3 

0.2168 

3 

4 

5.8238 

1 .0034 

3 

0.2174 

2 

4 

4.2443 

1.0069 

3 

S  -  6.0193 
DEADTIME  CORRECTION: 

M  AX.= 

-0.0116  MIN 

-0.0011  MAG. 


DIFFERENTIAL  MAGNITUDES  FOR  OUTPUT  FILE  870328.T04 


CYCLE-  l  FILTER=  4 

-2 

-3 

-4 

STAR  1- 

,  DAY=  28.1983  DMAG  = 

0.8486 

-0.7361 

-1.5069 

DAIR= 

-0.0004 

0.0001 

0.0001 

S  PAR  2- 

,  DAY=  28.1996  DMAG= 

-1.5846 

-2.3555 

DAIR= 

0.0006 

0.0006 

STAR 3- 

,  DAY=  28.2009  DMAG  = 

-0.7708 

DAIR= 

0.0000 

CYCLE=  2  FILTER=  2 

-2 

-3 

-4 

STAR  1- 

,  DAY=  28.2037  DMAG= 

0.7751 

-0.5203 

-1.3695 

DAIR= 

-0.0008 

0.0005 

0.0004 

STAR  2- 

,  DAY=  28.2049  DMAG  = 

-1.2955 

-2.1446 

DAIR= 

0.0012 

0.0011 

STAR  3- 

,  DAY=  28.2061  DMAG  = 

-0.8491 

DAIR= 

-0.0001 

-  100  - 


CYCLE=  3  FILTER=  2 

-2 

-3 

-4 

STAR  1- 

,  DAY=  28.2087  DMAG  = 

0.7698 

-0.5248 

-1.3730 

DAIR= 

-0.0010 

0.0008 

0.0006 

STAR  2- 

,  DAY=  28.2099  DMAG  = 

-1.2946 

-2.1428 

DAIR= 

0.0018 

0.0016 

STAR  3- 

,  DAY=  28.2112  DMAC«  = 

-0.8482 

DAIR  = 

-0.0002 

CYCLE=  4  FILTER=  4 

-2 

-3 

-4 

STAR  1- 

,  DAY=  28.2140  D\1AG= 

0.8470 

-0.7336 

-1.5093 

DAIR= 

-0.0014 

0.0011 

0.0008 

STAR  2- 

,  DAY=  28.2152  DMAG  = 

-1.5806 

-2.3563 

DAIR= 

0.0024 

0.0022 

STAR  3- 

.  DAY=  28.2165  DMAG  = 

-0.7757 

DAIR= 

-0.0002 

AVERAGED  DIFFERENTIAL  MAGNITUDES 

FILTER=  4 

-2 

-3 

-4 

STAR  1- 

,  DAY=  28.2061  DMAG  = 

0.8478 

-0.7348 

-1.5081 

STAR 2- 

,  DAY=  28.2074  DMAG  = 

-1.5826 

-2.3559 

STAR 3- 

.  DAY=  28.2087  DMAG  = 

-0.7733 

AVERAGED  DIFFERENTIAL  MAGNITUDES 

FILTER=  2 

-2 

-3 

~4 

STAR  1- 

,  DAY  =  28.2062  DMAG= 

0.7725 

-0.5226 

-1.3712 

STAR  2- 

,  DAY=  28.2074  DMAG= 

-1.2950 

-2.1437 

STAR  3- 

,  DAY=  28.2086  DMAG  = 

-0.8487 

MIDTIME  OF  OBSERVATIONS  WAS  1987  3  28.20746 


